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PREFACE 


An  understanding  of  the  statistical  nature  of  aircraft  wake 
vortex  decay  can  lead  to  improved  airport  capacity.  This  report 
presents  statistical  data  on  the  decay  of  wake  vortex  strength 
measured  with  the  Monostatic  Acoustic  Vortex  Sensing  System  (MAVSS). 
The  data  might  be  used  to  refine  the  wake  vortex  aircraft  cate¬ 
gories  (Heavy,  Large,  and  Small),  which  presently  are  predicated  on 
gross  certificated  takeoff  weight  rather  than  directly  on  vortex 
behavior. 

The  authors  would  like  to  thank  Myron  Clark  for  reviewing 
early  drafts  of  this  report  and  for  his  many  comments  toward 
improving  the  clarity  of  an  inherently  complex  analysis. 
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1.  INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  operational  hazard  to  air  traffic  posed  by  aircraft  wake 
vortices  is  ultimately  limited  by  the  decay  of  the  strength  of  the 
organized  vortical  motion.  Experiments  have  shown  that  the  decay 
can  take  several  forms,  generally  occurs  rapidly  once  it  begins, 
and  is  a  consequence  of  an  instability  in  the  vortex  configuration. 

Persistence  of  the  wake  vortex  hazard  depends  upon  aircraft 
parameters  (wingspan,  weight,  configuration,  engine  location,  etc.), 
meteorological  parameters  (wind  velocity,  wind  shear,  turbulence, 
atmospheric  stability,  pressure,  etc.),  and  decay  processes 
(vortex  linking,  bursting,  and  turbulent  diffusion).  Since  the 
decay  processes  occur  at  random  even  when  all  the  parameters  are 
fixed,  the  persistence  of  a  vortex  can  be  defined  only  through 
a  probability. 

The  current  designation  of  wake  vortex  separation  categories 
assigns  the  wake  vortex  hazard  to  a  single  aircraft  parameter, 
the  maximum  certificated  gross  takeoff  weight.  Of  necessity, 
this  simplified  procedure  gives  only  a  rough  indication  of  the 
wake  vortex  hazard.  The  actual  hazard  persistence  for  a  specified 
aircraft  has  a  wide  spread  because  of  variations  in  the  aircraft 
parameters,  variation  in  the  meteorological  conditions,  and  the 
stochastic  nature  of  vortex  decay. 

It  is  envisioned  that  eventually  aircraft  categories  will  be 
set  based  on  an  understanding  of  what  aircraft  parameters  in 
addition  to  weight  should  be  included.  Wingspan  must  certainly 
play  a  role,  and  engine  placement  is  also  likely  to  be  important. 
This  report  examines  these  issues.  If  vortex  decay  can  be 
characterized  by  a  few  aircraft  parameters,  then  one  can  hope 
to  have  a  more  rational  system  of  dealing  with  the  wake  vortex 
hazard. 
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Between  July  1976  and  September  1977  TSC  collected  data 
at  Chicago's  O'Hare  International  Airport  on  the  decay  of  wake 
vortex  strength.  The  Monostatic  Acoustic  Vortex  Sensing  System 
(MAVSS) ,  the  remote  sensing  technique  employed,  makes  use  of 
sound  energy  scattered  from  temperature  fluctuations  to  measure 
the  velocity  profile  of  a  vortex.  Volume  I  of  this  report  (Ref. 

1)  described  the  hardware  and  data  processing  involved  in  these 
measurements.  The  data  were  collected  during  normal  airport 
landing  operations.  Volume  II  (Ref.  2)  described  the  analysis  of 
whether  landing  B-707  and  DC-8  aircraft  need  to  be  divided  into 
Heavy  and  Large  categories  on  the  basis  of  their  wake-vortex 
hazard.  The  results  of  the  study  indicate  that  all  landing  B-707 
and  DC-8  aircraft  may  be  included  in  the  Large  wake-vortex  category. 
Volume  III  (Ref.  3)  summarized  the  results  of  Volume  II  in  terms 
of  safety  implications  of  categorizing  all  landing  B-707s  and 
DC-8s  as  Large  aircraft.  This  report  (Volume  IV  in  the  four- 
volume  series)  describes  the  various  statistical  methods  used  to 
understand  wake  vortex  decay  and  presents  data  on  all  common  jet 
transport  aircraft. 

Section  2  outlines  the  data  processing  procedures  and  de¬ 
scribes  in  detail  the  methods  used  to  analyze  the  MAVSS  data  and 
their  limitations.  Section  3  introduces  statistical  models  of 
vortex  decay.  Section  4  discusses  a  hazard  model  which  is  used 
to  relate  the  MAVSS  vortex  strength  measurements  to  an  encounter 
hazard  for  a  following  aircraft.  Section  5  relates  the  wake  vor¬ 
tex  characteristics  of  an  aircraft  to  its  physical  parameters 
such  as  wingspan,  weight,  engine  placement,  etc.  Empirical  models 
of  vortex  decay  are  developed;  two  theoretical  models  extensively 
used  in  the  literature  are  shown  to  be  unrealistic.  In  Section  6 
a  preliminary  look  at  recategorizing  aircraft  based  on  vortex 
behavior  is  presented. 

1 . 2  SUMMARY 

The  maximum  certificated  gross  weight  is  not  a  particularly 
good  parameter  for  classifying  aircraft  as  vortex  generators  or 
vortex  encounterers .  Wingspan  is  the  best  parameter  for 


1-2 


characterizing  a  vortex  encounterer.  Wingspan  with  some  correc¬ 
tions  for  engine  placement  can  also  be  used  to  characterize  jet 
transports  as  vortex  generators.  A  combination  of  weight  and 
wingspan  is  likely  to  be  needed  to  characterize  vortex- generating 
aircraft  with  wing  loading  significantly  different  from  that  of 
jet  transports. 

The  decay  of  wake  vortex  strength  can  be  represented  by  a 
simple  model  which  contains  a  few  aircraft-dependent  parameters. 

The  decay-time  parameter  is  observed  to  depend  much  less  upon 
aircraft  size  than  has  normally  been  assumed. 

The  modeled  vortex  decay  was  used  to  evaluate  the  safe 
separation  vortexwise  between  pairs  of  aircraft  under  the  assump¬ 
tion  that  the  hazard  probabilities  for  the  current  separation 
standards  are  acceptable  for  aircraft  pairs  which  occur  frequently. 
The  resulting  safe  separations  are  then  used  to  evaluate  briefly 
a  four-category  separation  system.  This  system  appears  to  offer 
the  possibility  of  improved  safety  while  simultaneously  relieving 
congestion  at  the  major  hub  airports.  Recommendations  for  further 
analysis  are  presented,  if  the  four-category  system  is  judged  to 
have  operational  merit. 
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2.  MAVSS  MEASUREMENTS 


The  method  of  processing  the  Monostatic  Acoustic  Vortex 
Sensing  System  (MAVSS)  data  is  described  in  detail  in  Volume  I 
of  this  report  (Ref.  1).  Only  a  summary  will  be  presented  here. 

The  methods  used  to  analyze  the  MAVSS  data  and  their  limitations, 
however,  will  be  examined  in  detail  herein. 

A  single  MAVSS  antenna  as  used  here  is  equivalent  to  a  60 
(or  90-)  meter  tower  with  16  evenly-spaced  vertically-pointing 
single-axis  anemometers.  It  measures  the  vertical  component  of 
the  wind  at  16  levels  or  range  gates  every  0.4  (or  0.6)  seconds 
with  a  vertical  resolution  of  3.5  (or  5.2)  meters  and  a  horizontal 
resolution  of  about  2  meters.  The  90-meter  range  (0.6-second 
sample  rate)  was  used  for  aircraft  taking  off  (Ref.  4).  The 
remote  sensing  capability  of  the  MAVSS  provides  two  major  advantages 
over  the  equivalent  instrumented  tower:  1)  it  does  not  protrude 
into  the  airspace  so  as  to  pose  a  safety  hazard  to  aircraft,  and 
2)  it  does  not  influence  the  decay  of  the  vortices  being  measured. 

The  actual  MAVSS  installation  consists  of  an  array  of  up 
to  10  antennas  positioned  along  a  baseline  perpendicular  to  the 
aircraft  flight  path.  A  measurement  of  vortex  strength  is  made 
whenever  a  vortex  drifts  over  an  antenna  in  the  array.  When  the 
ambient  wind  is  appropriate,  a  vortex  will  drift  past  many  antennas 
and  the  reasurements  will  show  how  the  vortex  decays  with  time. 
Figure  la  shows  the  MAVSS  trajectories  of  the  two  counterrotating 
vortices  generated  by  an  aircraft  taking  off.  The  vertical  ve¬ 
locity  signatures  of  the  vortices  are  shown  in  Figure  lb  for  each 
vortex  detection  at  an  antenna.  The  vortex  height  is  determined 
by  the  range  gate  showing  the  highest  velocities.  The  vortex 
arrival  time  is  indicated  by  the  change  in  sign  of  the  velocity. 

The  transport  velocity  of  a  vortex,  measured  by  the  arrival  times 
at  successive  antennas,  is  used  to  convert  the  velocity  time 
histories  to  the  spatial  velocity  profiles  in  Figure  lb.  One  of 
the  primary  reasons  for  the  success  of  the  MAVSS  is  that  it  is 
insensitive  to  the  ambient  wind  which  is  horizontal  near  the 
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FIGURE  1.  MAVSS  DATA:  (a)  VORTEX  TRAJECTORIES,  (b)  VELOCITY  SIGNATURES 


ground.  Any  departure  of  the  vertical  velocity  from  zero  can  be 
attributed  to  the  aircraft  wake.  Moreover,  the  vertical  velocity 
profile  at  the  vortex  height  is  a  good  measurement  of  the  vortex 
tangential  velocity  profile  under  the  assumption  of  axial  symmetry. 

Appendix  A  studies  the  detection  threshold  of  the  MAVSS  which 
plays  an  important  role  in  the  analysis  of  vortex  decay.  Ref¬ 
erence  1  contains  background  information  on  the  MAVSS  and  Appendix 
C  of  Reference  5  addresses  the  error  limitations  of  the  MAVSS. 

2.1  VORTEX  STRENGTH 

The  velocity  profile  V(r),  where  r  is  the  vortex  radius, 
measured  by  the  MAVSS  can  be  converted  to  the  circulation  or 
strength  profile  F(r)  by  the  relationship: 

T (r)  =  2irr  V(r)  .  (1) 

Figure  2  shows  the  relationship  between  V(r)  and  F(r)  for  an  ideal 
vortex  which  is  described  by  the  following  circulation  profile: 

F  (r)  =  TJ(1  +  Oc/r)2),  (2) 

where  is  the  circulation  at  large  radii  and  rc  is  the  core 
radius  which  represents  the  point  where  V(r)  is  maximum  and  also 
where  r(r_)  is  half  of  r  . 

v.  00 

The  parameter  used  to  evaluate  the  vortex  hazard  is  the 
average  circulation, 

f (r)  =  I  !  T  (r ' )  dr’.  (3J 

•'o 

The  integral  can  be  carried  out  for  the  form  in  Equation  2  to 
y  ie  Id 


r'(r)  =  r^ll  -  (rc/r)  tan' 1 (r/rc) ] . 


(4) 
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Figure  2  also  shows  how  r'(r)  compares  to  V(r)  and  T(r). 
r*(r)  approaches  the  limiting  value  F^  much  more  slowly  than 
does  r(r).  The  average  value  is  kept  below  by  the  region  of 
low  T(r)  below  about  2rc. 

The  MAVSS  actually  makes  two  measurements  of  V(r),  one  on 
either  side  of  the  vortex  core  (see  Figure  lb).  To  eliminate  some 
systematic  errors,  values  of  F’  are  accepted  only  if  the  data  from 
both  sides  can  be  averaged.  The  values  of  F*  for  four  values  of 
averaging  radius  (5,  10,  20,  and  30  meters)  are  stored  in  a  data 
base.  Values  for  15-meter  averaging  radius  are  subsequently 
obtained  by  interpolation. 

The  use  of  average  circulation  to  characterize  vortex  strength 
has  the  advantage  of  giving  a  more  stable  value  compared  to  other 
ways  of  processing  the  data,  such  as  making  a  least-squares  fit 
to  a  model  such  as  Equation  2  (Ref.  13).  The  stored  values  of  T' 
can  be  used,  nevertheless,  to  estimate  effective  values  of  T  and 
rc<  This  process  is  particularly  useful  if  T'(r)  is  averaged 
over  many  cases  before  and  rc  are  estimated.  If  F'(r)  is  known 
for  two  values  of  r  (a  and  b) ,  then  one  can  calculate  the  average 
of  the  circulation  between  radius  a  and  radius  b: 


r,(a,b)  =  jfb  T(r)  dr 

=  (br'(b)  -  ar'(a))/(b-a). 


(5) 


If  the  core  radius  is  5  m  or  less,  then  the  integral  between 
10  and  20  m  gives  a  reasonable  estimate  of  Ym.  This  estimate  can 
be  improved  if  the  core  radius  is  estimated.  Equation  4  can  be 
used  to  evaluate  the  ratio  of  the  two  average  circulations: 


r' 

-p- 


(a) 

m 


1  -  (rc/a)  tan_1(a/rc) 
1  -  (rc/b)  tan”*(b/rc) 


(6) 
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This  equation  can  be  solved  graphically  for  rc  (see  Figure  3) . 
One  can  then  use  rc  in  Equation  4  to  estimate  r^.  This  method 
of  analyzing  the  data  depends  primarily  on  the  area  of  the  circula¬ 
tion  deficit  for  r  <  2rc  and  is  not  sensitive  to  the  actual  shape 
of  the  vortex  circulation  profile. 


This  method  of  using  rc  to  represent  the  area  of  the  r(r) 
deficit  for  r  <  2rc  can  also  be  used  to  correct  for  one  of  the 
MAVSS  systematic  errors,  namely  an  underestimate  of  the  vortex 
velocity  for  small  vortex  radii.  Four  effects  contribute  to 
this  error:  (1)  The  coarse  spatial  resolution  of  the  MAVSS 
smears  out  the  center  of  the  vortex.  (2)  The  very  high  velocities 
(up  to  60  m/sec)  present  in  some  vortices  cannot  be  detected 
because  the  MAVSS  bandwidth  and  processing  cut  off  any  Doppler 
signals  above  approximately  13  m/sec.  (3)  The  laminar  flow  in 
vortex  cores  tends  to  suppress  the  turbulence  which  is  responsible 
for  producing  the  MAVSS  signals.  (4)  For  rapidly  moving  vortices 
the  sampling  rate  is  too  slow  (0.4  or  0.6  seconds/sample)  to 
accurately  characterize  the  vortex  core;  for  10  m/sec,  the  highest 
transport  velocity  accepted  for  vortex  analysis,  the  sample  spac¬ 
ing  is  4  or  6  meters.  The  loss  in  MAVSS  signal  at  small  radii 
will  increase  the  area  of  the  T(r)  deficit  and  lead  to  a  larger  rc 
but  with  the  same  value  of  r  .  If  r_„  is  the  actual  core  radius, 
then  the  actual  average  circulation  r'(r)  is  related  to  the  measured 

a 

values  of  F'(r)  and  rc  by 


r’(r) 


r'(r) 


d  -  (rca/r)  tan_1(r/rca)) 
(1  -  (rc/r)  tan"1(r/rc)) 


(7) 


Table  1  illustrates  the  sorts  of  errors  which  can  arise  for 
some  estimates  of  core  radius  to  be  encountered  in  subsequent 
analysis.  The  corrections  can  be  substantial  especially  for  5- 
meter  average  radius.  The  measured  values  of  rc  (shown  later 
in  Table  10)  are  roughly  5.5  m  for  widebody  aircraft  and  4.0  m 
for  other  aircraft  types.  The  core  size  for  small  aircraft  with 
no  wing-mounted  engines,  such  as  the  DC-9  and  the  B-727,  could 
be  as  small  as  0.5  meters.  Aircraft  with  four  wing-mounted  engines 
tend  to  have  larger  core  sizes. 
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FIGURE  3.  DEPENDENCE  OF  r,(a)/r,(b)  UPON  rr/a  FOR  TWO  VALUES  OF  a/b 
THE  SCALE  CHANGES  FROM  LINEAR  TO  LOGARITHMIC  AT  r  /a  =  0.2. 


TABLE  1.  AVERAGE  CIRCULATION  CORRECTIONS  FOR  ERRORS  IN  CORE 
RADIUS 


r^(r)/r»(r) 


;a  (m) 

rc  (m) 

r  =  5  m 

r  =  10  m 

r  e  15  m 

r  =  20  m 

0.5 

4.0 

3.01 

1.76 

1.33 

1.5 

4.0 

2.18 

1.50 

1.22 

2.5 

4.0 

1.58 

1.28 

1.13 

2.5 

5.0 

2.37 

1.62 

•1.38 

1.28 

• 

2.2  VORTEX  STRENGTH  HISTORY 

The  primary  method  selected  for  analyzing  the  MAVSS  data  re¬ 
quires  that  the  time  history  of  the  vortex  strength  (i.e.,  average 
circulation)  be  known  for  each  vortex  measured.  The  MAVSS  data 
base  contains  strength  values  only  at  times  when  a  vortex  passes 
over  an  antenna  in  the  array.  These  values  are  processed  by  a 
combination  of  interpolation  and  extrapolation  techniques  to  obtain 
a  time  history  of  the  vortex  strength.  The  following  assumptions 
are  made: 

1.  The  vortex  strength  is  constant  until  the  vortex  is 
first  detected. 

2.  The  vortex  strength  is  obtained  by  interpolation  between 
successive  vortex  detections. 

3.  The  vortex  strength  is  assumed  to  be  zero  if  it  cannot  be 
detected  at  the  next  antenna  after  the  last  detection. 

The  vortex  strength  is  then  extrapolated  to  zero  at 

the  time  it  would  be  expected  to  reach  the  next  antenna, 
based  on  the  latest  value  for  transport  velocity.  This 
extrapolation  is  excluded  if  the  last  detection  is  in  the 
last  antenna  in  the  array  or  if  the  expected  arrival  time 
is  so  close  to  the  next  aircraft  landing  that  the  next 
detection  could  be  obscured  by  noise. 
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The  death  of  a  vortex  in  3)  is  signaled  by  not  being  detected  in 
an  antenna  over  which  it  passes.  The  two  exclusions  on  the 
extrapolation  in  3)  are  designed  to  eliminate  any  causes  for  non¬ 
detection  (e.g.,  aircraft  noise)  other  than  decay  to  a  strength 
below  the  detection  threshold.  Appendix  A  evaluates  the  detec¬ 
tion  threshold.  Section  2.1  of  Volume  II  (Ref.  2)  illustrates 
this  vortex  history  procedure  for  actual  data.  Table  2  lists  the 
number  of  vortex  deaths  measured  for  each  aircraft  type. 


TABLE  2.  NUMBER  OF  MEASURED  VORTEX  DEATHS 


Aircraft  Type 

Vortex  1 

Vortex  2 

DC-9 

318 

171 

B-  7  37 

89 

50 

B-727 

848 

464 

B-707 

249 

162 

DC-3 

133 

87 

DC-8H 

75 

51 

B-707H 

31 

17 

DC-10 

261 

194 

L- 1011 

62 

39 

B-747 

135 

69 

2.3  VORTEX  STRENGTH  ERRORS 

In  the  vortex  history  the  two  extrapolations  1)  and  3)  are 
more  susceptible  to  error.  In  the  case  of  a  slowly  moving  vortex 
the  first  detection  can  occur  after  a  significant  amount  of  decay. 
In  recent  data  analysis,  assumption  1)  has  been  modified  to  allow 
extrapolation  to  zero  time  only  if  the  first  detection  occurs 
before  a  vortex  age  of  20  seconds.  The  assignment  of  zero  strength 
to  a  missed  detection  in  3)  is  also  a  rather  crude  approximation 
in  view  of  the  finite  detection  threshold  and  the  total  uncertainty 
in  the  actual  death  position  of  a  vortex.  The  errors  associated 
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with  these  extrapolations  are  mitigated  by  the  requirement  that 
all  vortices  included  in  the  data  analysis  are  detected  in  at 
least  two  antennas.  This  requirement  also  assures  reasonable 
accuracy  since  the  vortex  strength  is  proportional  to  the 
measured  transport  velocity  which  can  be  reliably  estimated  from 
two  detections. 

An  other  possible  source  of  strength  error  can  enter  into 
the  vortex  strength  at  early  times  fe.g.,  less  than  20  seconds) 
when  the  vortex  is  still  descending  toward  the  ground.  The  use 
of  a  single  range  gate  for  the  tangential  velocity  profile  can 
underestimate  the  vortex  strength  when  the  maximum  velocities 
occur  at  different  heights  before  and  after  the  vortex  passes. 
Readings  as  much  as  25  percent  low  have  been  seen  from  this  effect. 
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3.  STATISTICS  OF  VORTEX  DECAY 


The  decay  of  wake  vortices  must  be  described  as  a  statistical 
process  for  two  reasons.  First,  the  condition  of  the  atmosphere 
depends  upon  too  many  variables,  (e.g.,  wind  distribution,  tempera¬ 
ture  distribution,  etc.)  to  be  specified  exactly.  Since  the 
atmospheric  conditions  have  a  strong  influence  on  vortex  decay, 
a  statistical  description  is  needed  to  account  for  unknown  vari¬ 
ables.  Second,  many  vortex  decay  mechanisms  are  themselves  random 
processes  which  cannot  be  defined  by  the  average  properties  of 
the  atmosphere.  A  statistical  description  of  vortex  decay  is 
therefore  the  best  that  can  be  expected. 

Unfortunately,  the  average  decay  time  of  wake  vortices  is  of 
little  use  in  assessing  the  wake  vortex  hazard.  Assured  safe 
separation  for  a  following  aircraft  requires  a  determination  of 
the  maximum  duration  of  the  wake  vortex  hazard.  Such  a  determina¬ 
tion  immediately  leads  to  the  difficult  problem  of  measuring  small 
probabilities . 

3.1  VORTEX  HAZARD  PROBABILITY 

Section  4.1  describes  a  method  for  defining  a  vortex  strength 
hazard  threshold  r|.  For  the  present  discussion  it  is  sufficient 
to  note  that  a  wake  vortex  poses  a  potential  hazard  to  a  following 
aircraft  separated  by  a  time  t  if  the  vortex  strength  r'(t)  at  age 
t  is  greater  than  Tj  and  poses  no  hazard  if  r'(t)  is  less  than  Tj. 
The  vortex  hazard  probability  F(r^,t)  predicts  what  fraction  of 
the  wake  vortices  generated  will  have  a  strength  less  than  r^.  at 
age  t.  The  dependence  of  F(F-J,,t)  upon  strength  threshold  and  age 
characterizes  the  decay  of  the  wake  vortex  hazard. 

Probabilities  can  be  measured  in  a  simple  way.  The  probabil¬ 
ity  of  a  certain  event  is  given  by  the  ratio  of  the  number  of  times 
the  event  occurs  divided  by  the  total  number  of  times  it  could 
have  occurred.  Because  of  statistical  fluctuations,  an  accurate 
measurement  of  probability  requires  a  large  number  of  events, 
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particularly  if  the  probability  is  small.  For  example,  if  IT 
events  are  expected  in  1000  cases,  the  probability  F  is  IT/1000. 

For  IT  much  less  than  1000,  the  root-mean-square  (rms)  variation 
in  number  of  events  for  different  selections  of  1000  cases  is 
FT/ Thus,  the  rms  variation  in  F  becomes  large  for  IT  less  than 
10.  In  the  particular  case  of  F(r|,t),  the  vortex  hazard 
probability  at  a  time  t,  F  is  given  by  the  ratio  of  the  number  of 
cases  with  strength  equal  to  or  greater  than  T.J.  at  time  t  divided 
by  the  total  number  of  cases  with  valid  strength  measurements  at 
t ime  t . 

3.2  HAZARD  PROBABILITY  DECAY 

Figure  4  shows  an  example  of  how  the  vortex  hazard  probabil¬ 
ity  decays  in  time.  The  vortex  time  histories  (Section  2.2)  were 
evaluated  for  all  B-707  vortices  with  a  10-meter  averaging  radius 
and  a  hazard  threshold  of  100  m  /sec.  The  probability  was 
calculated  at  10-second  intervals.  The  logarithm  of  the  hazard 
probability  is  observed  to  decay  more  rapidly  as  the  vortex  age 
increases.  In  this  plot  the  total  number  of  vortices  measured  is 
about  500.  A  probability  of  0.002  corresponds  to  one  case  above 
hazard  threshold  while  a  probability  of  0.02  corresponds  to  10 
cases  above  the  threshold.  The  earlier  discussion  of  statistical 
fluctuations  indicates  that  the  measured  probabilities  below  0.02 
can  have  large  statistical  errors  (see  Appendix  D  of  Volume  II, 
Ref.  2). 

Early  in  the  analysis  of  the  hazard  decay  plots  it  was  found 

that  the  decay  in  hazard  probability  F  becomes  approximately 

2 

linear  when  plotted  as  log  F  vs.  t  .  Figure  5  shows  F(T2.,t)  for 

7  1 

strength  thresholds  of  150,  100,  75,  50,  and  30  m  /sec  plotted 
2 

against  t  .  The  second  curve  from  the  left  is  the  data  of  Figure 
4.  The  probability  decreases  as  the  threshold  increases.  The 
format  of  Figure  5  has  become  the  standard  display  used  to  analyze 
the  vortex  decay  data.  It  shows  in  a  single  picture  how  the 
vortex  hazard  probability  depends  upon  the  vortex  age  and  the 
assumed  hazard  threshold.  Six  hazard  strength  values  (the  five 
above  and  200  m  /sec)  were  used  in  the  standard  display  of  vortex 
decay . 
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FIGURE  4.  HAZARD  PROBABILITY  VERSUS  VORTEX  AGE  FOR  ALL  B-707 
VORTICES  FOR  T*  =  100  mZ/sec  AND  10-m  SEMISPAN 
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Some  of  the  details  of  Figure  5  require  explanation.  The 
x's  plotted  in  Figure  5  correspond  to  the  probability  correspond¬ 
ing  to  one  case.  Since  zero  cannot  be  plotted  on  a  logarithmic 
scale,  the  plots  are  terminated  at  the  level  of  one  case.  If  the 
probability  drops  to  zero  from  a  value  higher  than  one  case,  a 
dashed  line  is  drawn  to  the  one-case  value  at  the  age  where  zero 
occurred.  The  age  intervals  are  10  seconds  as  in  Figure  4.  The 
single-case  probability,  marked  by  x's  in  Figure  5,  rises  with 
vortex  age  because  the  number  of  valid  strength  histories  de¬ 
creases  with  vortex  age.  Vortices  are  lost  by  passing  the  end 
of  the  antenna  array  or  by  interference  from  the  noise  of  the 
following  aircraft.  The  number  of  cases  reaches  a  steady  level 
after  about  80  seconds,  when  most  of  the  valid  measurements  are 
those  of  vortices  which  have  been  declared  dead  according  to  the 
restrictions  of  Section  2.2. 

A  straight-line  fit  to  the  curves  of  Figure  5  corresponds 
to  the  functional  form 


F(r^,t)=  exp  (-a(t2-d))  ,  (8) 

2 

where  a  represents  the  slope  and  d  the  delay  in  t  to  the  point 
where  F  =  1.  This  functional  dependence  describes  the  vortex 
decay  for  more  than  three  decades  for  the  aircraft  with  the  most 
data  (notably,  the  B-727).  This  form  also  appears  to  give  a 
reasonable  asymptotic  fit  to  the  data  no  matter  how  the  data  are 
disaggregated  according  to  meteorological  conditions. 

The  functional  dependence  of  equation  8  implies  a  decay 
probability  that  increases  in  time,  which  is  plausible  for  the 
various  modes  of  vortex  decay.  The  probability  of  vortex  linking 
and  core  bursting  is  observed  to  increase  with  time  (Refs.  6  and 
7).  The  vortex  decay  due  to  turbulent  diffusion  may  also  in¬ 
crease  with  time  as  vortex  energy  is  converted  to  turbulent  energy 
and  as  vorticity  is  diffused  away  from  the  vortex  core  to  regions 
where  it  can  annihilate  with  vorticity  of  opposite  sign. 
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The  fact  that  the  vortex  hazard  decays  as  exp(-at  )  accounts 
for  the  observed  safety  of  the  air  transportation  system.  If  the 
decay  were  slower,  for  example,  as  expC-t/x),  the  accident  rate 
could  not  remain  small  except  with  inordinately  large  aircraft 
separations.  Suppose  the  hazard  probability  decayed  exponentially 
with  a  time  constant  t  as  short  as  30  seconds  (which  is  much 
shorter  than  the  e decay  for  many  of  the  aircraft  at  reasonable 
hazard  strength  levels;  it  would  require  138  seconds  for  the 
hazard  probability  to  drop  to  1  percent  and  207  seconds  to  drop 
to  0.1  percent. 


3.3  VORTEX  STRENGTH  DISTRIBUTION 

This  section  develops  the  mathematical  formalism  appropriate 
to  a  probabal istic  description  of  vortex  decay  and  then  examines 
how  the  vortex  strength  distribution  varies  in  time.  Readers 
interested  only  in  hazard  decay  should  move  ahead  to  Section  4.1. 

At  any  given  time  the  vortex  strength  T '  can  be  described  by 
the  distribution  function  P(r',t)  which  gives  the  probability 
P  df  of  the  vortex  strength  lying  in  a  small  interval  dF'.  The 
normalization  of  P  is  given  by: 


/ 


p (r *  ,t)  dr’ 


0 


1. 


(9) 


A  knowledge  of  P(r',t)  constitutes  a  complete  description  of 
vortex  decay  under  the  conditions  specified  for  the  distribution 
function  (e.g.,  meteorological  conditions,  aircraft  type,  air¬ 
craft  weights,  etc.).  For  example,  the  average  vortex  strength 
r * ,  is  given  by 


r'  (t) 


f* P(r\t)  r’ 

o 


dr* . 


(10) 


I 


* 


3-6 


The  probability  F(r|,t)  that  the  vortex  strength  is  above  the 
hazard  threshold  is  given  by 

F(r{,t)  =  fm  P(r*  ,t)  dr*.  (11) 

rj. 

1 

The  evaluation  of  the  vortex  strength  distribution  function 
P(T’,t)  is  similar  to  that  of  the  probability  F(r^.,t).  The  vortex 
history  algorithms  of  Section  2.2  are  used  to  get  the  strength  of 
each  vortex  at  age  t.  The  vortices  are  then  sorted  into  bins  20 
m  /sec  wide.  The  strength  distribution  at  the  midpoint  of  the 
bin  is  then  defined  as  the  ratio  of  the  number  of  vortices  in  the 
bin  divided  by  the  bin  width  and  the  total  number  of  valid  strength 
measurements  at  the  time  t.  Figure  6  shows  the  resulting  strength 
distribution  for  all  B-707  landing  vortices  at  a  number  of  different 
ages.  The  initial  distribution  (10  sec)  consists  of  a  well- 
defined  (40  percent  full  width  at  half  maximum)  peak  with  a  tail 
extending  to  low  strengths.  Presumably  the  tail  is  due  to 
vortices  which  have  decayed  before  they  were  first  detected.  The 
peak  in  the  distribution  function  P  corresponds  to  the  most  probable 
vortex  strength.  The  initial  most  probable  strengths  are  listed 
in  Table  3  for  the  major  aircraft  types  using  all  the  data  col¬ 
lected. 

The  vortex  decay  in  Figure  6  results  in  a  steady  decrease 
in  the  most  probable  vortex  strength  and  an  increase  in  the  width 
of  the  distribution  function.  Eventually,  a  large  portion  of  the 
distribution  function  drops  below  the  detection  threshold  of  30 
or  40  m  /sec  (Appendix  A).  Figure  7  shows  the  probability  of  the 
vortex  strength  remaining  nonzero  as  a  function  of  time.  The  long 
duration  of  nonzero  strength  compared  to  Figure  4,  for  example, 
is  due  to  the  nature  of  the  extrapolation  process  by  which  a 
vortex  is  declared  defunct  (Section  2.2).  The  extrapolation 
procedure  leads  to  the  sizable  vortex  distribution  below  the 
detection  threshold  in  Figure  6. 


TABLE  3.  MOST  PROBABLE  INITIAL  VORTEX  STRENGTH  (m2/sec) 
FOR  LANDING  AIRCRAFT 


Aircraft 

Averaging 

Radius 

5  m 

10m 

15  m 

20  m 

DC-9 

46 

85 

95 

110 

B-737 

41 

85 

98 

x07 

B-727 

56 

105 

120 

140 

B-  707 

76 

132 

146 

165 

DC-8 

76 

135 

155 

185 

DC-8H 

72 

135 

145 

170 

B-707H 

80 

135 

150 

175 

DC-10 

60 

135 

165 

195 

L- 1011 

65 

140 

175 

205 

B-  747 

90 

185 

225 

260 

Consider  how  the  features  of  Figure  5  are  produced  by  the 
decay  in  the  strength  distribution  of  Figure  6.  The  initial 

2 

probability  is  near  one  when  the  threshold  is  below  the  135  m  /sec 
peak  in  the  initial  distribution.  It  is  much  reduced  when  the 

7 

threshold  is  above  the  peak,  as  for  150  m  /sec.  Another  feature 

of  the  decay  is  that  a  progressively  increasing  delay  in  the  decay 

occurs  when  the  threshold  is  dropped  farther  and  farther  below  the 

2 

initial  strength.  After  the  delay  the  log  F  versus  t  curve  again 
becomes  roughly  linear.  This  delay  represents  the  time  it  takes 
for  the  peak  in  the  distribution  curve  to  drop  below  the  thres¬ 
hold. 

A  number  of  different  methods  were  developed  to  model  the 
vortex  decay  of  Figure  6.  They  are  described  in  Appendix  B. 

Some  of  them  will  be  used  in  later  sections  of  this  report. 
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4.  hazard  analysis 


A  determination  of  the  duration  of  the  wake  vortex  hazard  is 
a  critical  element  in  evaluating  the  safety  of  wake  vortex  separa¬ 
tions.  A  hazard  model  is  needed  to  relate  the  MAVSS  strength 
measurements  of  Section  2  to  the  encounter  hazard  for  a  following 
aircraft.  A  number  of  improvements  to  the  traditional  hazard 
model  will  be  presented.  The  amount  of  MAVSS  data  collected  is 
sufficient  to  make  estimates  of  the  statistics  of  the  hazard 
duration  (as  outlined  in  Section  3).  In  addition,  the  dependence 
of  hazard  duration  upon  the  meteorological  conditions  will  be 
discussed. 

4.1  ROLLING  MOMENT  HAZARD  MODEL 

Volume  II  (Section  4.2.1)  of  this  report  (Ref.  2)  derives 
the  following  expression  for  the  vortex  strength  hazard  threshold: 

rf(b/2)  =  1  KfbVp,  (12) 

where  the  parameters  of  the  following  aircraft  are  b  =  the  wingspan, 
V  =  the  airspeed,  and  p  =  the  maximum  nondimens ional  roll  rate. 

This  hazard  model  assumes  that  a  vortex  poses  a  hazard  if  its 
induced  rolling  moment  is  greater  than  a  fraction  f  of  the  maximum 
roll  control  of  the  aircraft.  The  derivation  of  Equation  12  assumes 
that  the  value  of  P  is  independent  of  the  shape  of  the  vortex 
velocity  profile.  Equation  12  with  K  =  1  is  correct  for  a  velocity 
profile  which  increases  linearly  with  radius.  Appendix  B  of 
Reference  5  uses  vortex-lattice  theory  to  calculate  the  corrections 
for  other  velocity  profiles.  Table  4  lists  the  correction  factors 
K  for  two  aircraft  types  and  several  velocity  profiles:  1)  con¬ 
stant  velocity  and  2)  Equations  1  and  2  with  rc  =  0,  2.5,  and 
5  meters.  Deviations  from  K  =  1  are  larger  for  the  DC-9  because 
of  greater  taper  (and  perhaps  sweep)  in  the  wing  planform. 
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TABLE  4.  INDUCED  ROLLING  MOMENT  CORRECTION  FACTORS 


AIRCRAFT 

WINGSPAN 

CORRECTION  FACTO 

R  K 

(m) 

rc=0 

rc=2 . 5 

rc=5  m 

T-  37 

10.3 

.92 

.  82 

.91 

.95 

DC-9  (flaps) 

28.4 

.87 

.71 

.78 

.85 

The  basic  equation  (Equation  12)  of  the  hazard  model  can  be 
used  to  determine  the  wingspan  dependence  of  the  vortex  hazard. 
Dividing  Equation  12  by  b/2  isolates  the  wingspan  dependence: 

2l\j.(b/2)/b  =  Ht  =  KfVp.  (13) 

The  vortex  model  (Equation  4)  can  be  used  to  evaluate  r'(r)/r  as 
a  function  of  iOT  and  rc.  The  results  are  plotted  in  Figure  8. 

The  range  of  wingspans  experiencing  a  hazard  from  a  particular 
vortex  is  critically  dependent  on  the  value  of  HT  in  Equation  12. 

The  value  in  Figure  8  yields  no  hazardous  wingspans.  The 

value  Hj2  yields  a  hazard  for  wingspans  between  2  rc  and  5.2  rc . 
Reducing  H.^  by  a  factor  of  two  to  drastically  increases  the 

hazard  wingspan  range  to  .69  rc  through  14.8  r  .  For  further 
reductions  in  the  lower  hazard  bound  will  decrease  proportionally 
to  Hj  and  the  upper  bound  inversely  with  Hj. 

The  fact  that  a  given  vortex  is  safe  for  both  small  and  large 
wingspans  may  be  surprising  at  first  glance,  but  it  can  be  readily 
understood  by  careful  examination.  For  large  wingspans  the  vortex 
velocities  are  too  small  over  much  of  the  wingspan  to  present  a 
hazard.  For  wingspans  smaller  than  the  core  size  the  encountering 
aircraft  experiences  only  the  local  velocity  gradient  in  the 
core  and  not  the  full  velocity  variation  of  the  vortex. 


The  hazard  analysis  of  this  report  (also  Volume  II)  assumes 
a  hazard  threshold  of 


HT  =  10  f  (m2/sec) .  (14) 

This  value  is  consistent  with  the  parameters  K  =  1,  V  =  68  m/sec 
(130  knots),  and  p  =  0.07.  Table  5  shows  the  average  circulation 
hazard  threshold  values  from  Equation  14  with  f  =  1.0.  The 
sensitivity  of  the  results  to  the  parameter  f  will,  in  part, 
also  serve  to  illustrate  the  effects  of  varying  K,  V,  and  p. 
Variations  in  f  can  also  account  for  the  MAVSS  errors  estimated 
in  Table  1.  Possible  variations  in  K  are  illustrated  in  Table  3. 
Conservative  estimates  for  £  (Ref.  8)  are  0.06  for  jet  transports 
and  0.08  for  general  aviation  aircraft.  These  numbers  refer  to 
the  roll  control  of  ailerons  alone.  The  use  of  spoilers  for 
additional  roll  control  can  lead  to  p  =  0.12.  Military  aircraft 
generally  have  even  larger  values  of  p  (0.12  for  the  T-37  used  in 
vortex  encounter  flight  tests  and  0.18  for  the  F-86).  Because  of 
pilot  response  times,  the  full  roll-control  capability  of  an 
aircraft  cannot  be  used  to  compensate  for  a  vortex  encounter. 

Both  flight  encounters  and  simulator  studies  indicate  that  the 
actual  hazard  threshold  for  f  is  more  like  0.4  than  1.0. 

TABLE  5.  HAZARD  THRESHOLDS  FOR  f  =  1 


Aircraft  Type 

Semi  Span 

Strength  Threshold 

Small  GA 

5  m 

50  m2/sec 

Large  GA 

10  m 

100  m2/sec 

DC-9 

15  m 

150  m2/sec 

B-  707 

20  m 

200  m2/sec 

The  MAVSS  detection  thresholds  discussed  in  Appendix  A  lead 

to  a  limit  on  the  value  of  f  in  Equation 

14  which  can  be  examined 

experimentally.  The  approximate  detection  thresholds  for  5,  10, 
and  20-meter  averaging  radii  are,  respectively,  f  =  0.6,  0.4,  and 
0.3. 
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4.2  HAZARD  DECAY 


The  decay  of  the  vortex  hazard  probability  was  investigated 
in  Section  3.2  before  the  description  of  the  hazard  model  was  de¬ 
scribed  in  detail.  The  observations  there  can  now  be  interpreted 

in  terms  of  the  hazard  model  parameter  f.  The  choice  of  hazard 

2 

threshold  T-J.  =  100  m  /sec  for  10-m  averaging  radius  in  Figure  4 

corresponds  to  f  =  1.0  according  to  Equations  13  and  14.  In 

2 

Figure  5  the  curves  for  r^,  =  150,  75,  50,  and  30  m  /sec  correspond 
to  f  =  1.5  ,  0.75,  0.50  ,  and  0.30,  respectively.  The  different 
curves  thus  show  how  the  selected  value  of  f  affects  the  hazard 
decay.  As  one  would  expect,  the  hazard  lasts  longer  for  the  lower 
values  of  f. 

Figure  5  was  based  on  all  the  B-707  vortices  measured. 
Information  about  meteorological  effects  on  vortex  decay  can  be 
obtained  by  disaggregating  the  measurements  according  to  the 
ambient  wind  or  some  other  meteorological  parameter.  The  primary 
limit  on  disaggregating  the  data  is  the  necessity  of  retaining 
enough  cases  to  make  statistically  significant  measurements. 

The  most  significant  parameter  affecting  vortex  decay  was 
found  to  be  the  vortex  position  relative  to  the  ambient  crosswind. 
The  downwind  vortex,  which  is  the  first  to  reach  a  MAVSS  antenna 
(termed  vortex  1),  decays  more  rapidly  than  the  upwind  vortex,  which 
is  the  second  vortex  to  reach  an  antenna  (termed  vortex  2) . 

Figure  9  shows  the  data  of  Figure  5  disaggregated  into  vortex  1 
and  vortex  2.  The  operational  significance  of  this  effect  lies 
in  the  fact  that  vortex  2  is  the  one  which  could  linger  near  the 
extended  runway  centerline  and  pose  a  hazard  to  an  aircraft  fol¬ 
lowing  on  the  same  runway,  while  vortex  1  is  more  likely  to  drift 
toward  any  parallel  runway.  Most  of  the  subsequent  hazard 
analysis  (e.g.,  Section  6)  will  deal  with  the  decay  of  vortex  2 
which  determines  the  single  runway  vortex  separation  standards. 

Figure  10  shows  the  data  of  Figure  9  analyzed  with  the  recent 
algorithm  which  excludes  extrapolation  to  zero  time  unless  the 
first  vortex  detection  is  before  20  seconds.  The  result  of  this 
restriction  is  a  significant  reduction  in  the  number  of  cases  at 


VORTEX  HAZARD  PROBABILITY  (r* 


FIGURE  9.  COMPARISON  OF  VORTEX  DECAY  FOR  VORTEX  1  AND  VORTEX  2: 
B- 707  AIRCRAFT,  10-m  AVERAGING  RADIUS 


1 


early  times.  In  fact,  the  number  of  cases  becomes  relatively 
independent  of  vortex  age.  Appendix  C  presents  plots  like  Figure 
9  for  all  common  jet  transport  aircraft  and  averaging  radii. 

High  windspeeds  are  observed  to  increase  the  rate  of  vortex 
decay  (see  Appendix  C) .  The  functional  dependence  upon  wind  speed 
has  not  been  determined,  however,  because  the  effect  is  not  much 
greater  than  the  statistical  uncertainties. 


i 
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5.  AIRCRAFT  DEPENDENCE  OF  VORTEX  CHARACTERISTICS 


This  section  relates  the  wake  vortex  characteristics  of  an 
aircraft  to  its  physical  parameters  such  as  wingspan,  weight, 
engine  placement,  etc.  Vortex  data  are  available  for  American-made 
commercial  jet  transports.  Table  6  lists  the  relevant  parameters 
(Ref.  9)  for  various  models  of  these  aircraft. 

5.1  INITIAL  VORTEX  STRENGTH 

5 . 1 . 1  Theory 

The  classical  expression  for  the  dependence  of  total  wake 
vortex  strength  upon  aircraft  parameters  (Ref.  10)  is: 

=  cw/pbv,  (is) 

where  C  is  a  constant  depending  upon  the  wing  loading  distribution 
(C  -  4 / rr  for  elliptic  loading,  C  =  1  for  constant  loading),  W  is 
the  aircraft  weight,  p  is  the  air  density,  b  is  the  wingspan,  and 
V  is  the  airspeed.  The  total  circulation  increases  with  air¬ 
craft  weight  and  decreases  with  airspeed  and  wingspan.  Heavier 
aircraft  tend  to  have  larger  wingspans. 

The  simple  aircraft  dependence  of  Equation  15  is  in  the  three 
parameters  W,  b,  and  V.  The  wing-loading  distribution  factor  C 
is  much  more  difficult  to  address.  The  expression  W/bV  was  evaluated 
for  the  maximum  landing  weight  in  Table  6  for  all  the  aircraft 

models.  Two  values  were  used  for  V:  the  listed  landing  speed  and 
a  fixed  speed  of  150  mph.  Table  7  shows  the  values  for  W/bV 
assigned  to  each  type  as  a  mean  value  for  the  different  models  in 
use  at  the  time  of  our  data  collection.  A  mean  wingspan  is  also 
assigned.  The  spread  in  values  for  different  models  is  shown  in 
parentheses.  The  values  of  Table  7  are  converted  to  estimates  of 
FfX)  in  Table  8  through  the  use  of  Equation  15.  The  values  are  also 
converted  to  metric  units  in  order  to  allow  comparison  with  the 
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TABLE  6.  AIRCRAFT  MODEL  PARAMETERS 


MODEL 

WINGSPAN 

(ft) 

b 

BOEING 

707-120 

130.9 

707-120B 

130.9 

707-320 

142.4 

707-320B/C 

145.8 

707/420 

142.4 

720 

130.9 

720B 

130.9 

727-100 

108.0 

727-200 

108.0 

737-100 

93.0 

737-200 

93.0 

747-100B 

195.7 

747SR 

195.7 

747-200B/C 

195.7 

747SP 

195.7 

MCDONNEL  OO'JGLAS 


DC- 8- 10 

142.3 

DC-8-20 

142.3 

DC-8- 30/40 

142.3 

DC-8- 50 

142.3 

DC-8-61 

142.3 

DC-8-62 

148.4 

DC-8-63 

148.4 

DC-9- 10 

89.4 

DC-9- 20 

93.3 

DC-9- 30 

93.3 

DC-9-40 

93.3 

DC-9-50 

93.3 

DC-9-80 

107.8 

DC-10- 10 

155.3 

DC-10- 30 

165.3 

DC-10-40 

165.3 

LOCKHEAD-CAL I FORNIA 


TABLE  7.  AIRCRAFT  TYPE  PARAMETERS 


TYPE 

Wingspan  (Range) 

"(ft) 

WL/VLb  (Range) 

(lbs/  Mph,f  t; 

DC-9 

93.0  (89.4-93.3) 

8.3  (6. 3-8. 6) 

7.2  (6. 1-7. 9) 

B-737 

93.0 

8.3  (8. 1-8.4) 

7.3  (7. 2-7. 5) 

B-727 

108.0 

9.7  (9. 4-9. 9) 

9.3  (8. 8-9. 5) 

B-707 

130.9 

9.0  (8. 8-9. 2) 

9.3  (8. 9-9. 7) 

DC-8 

142.3 

9.2  (9. 2-9. 3) 

9.2  (9. 0-9. 4) 

DC-8H 

145  (142.3-148.4) 

10.4  (9.5-11.3) 

10.5  (9.7-11.2) 

B-707H 

144  (142.4-145.8) 

9.8  (9.0-10.7) 

10.5  (9.7-11.3) 

DC- 10 

160  (155.3-165.3) 

15.5  (15.1-15.8) 

15.9  (15.6-16.3) 

L-1011 

155.3 

14.3  (14.1-14.5) 

15.6  (15.4-15.8) 

B-747 

195.7 

17.6  (17.2-17.8) 

18.7  (17.9-19.2) 

I 
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TABLE  8.  INITIAL  VORTEX  STRENGTH 


AIRCRAFT 

TYPE 

WINGSPAN 

CALCULATED  1 

1 

Hiimi 

P (10-20  a) (4) 

mmm 

!H3» 

■ 

(m) 

KQZfiSHl 

(m^/sec) 

(m^/sec) 

(m^/sec) 

(m^/sec) 

1 

28.3 

270 

234 

179 

200 

161 

WSB^m 

28.3 

270 

237 

182 

203 

179 

B-727 

32.9 

315 

302 

223 

249 

230 

B-707 

39.9 

293 

302 

240 

259 

278 

DC-8 

43.4 

299 

299 

275 

298 

290 

DC-8H 

44.2 

338 

341 

268 

289 

281 

B-707H 

43.9 

319 

341 

238 

249 

286 

DC- 10 

48.8 

504 

517 

298 

348 

341 

L-1011 

47.3 

465 

507 

307 

354 

377 

B-747 

59.6 

572 

608 

417 

497 

499 

(1)  At  Maximum  Landing  Weight,  Temperature  *  0°C,  Altitude  *  sea  level, 
and  Elliptic  Wing  Loading. 

(2)  Published  Landing  Speed. 

(3)  Constant  Landing  Speed  (150  Mph) . 

(4)  Using  Vortex  Detections  Between  10  and  15  seconds. 

(5)  Linear  Fit  to  Vortex  Decay. 


I 


experimental  measurements.  The  air  density  p  was  taken  as  1.28 
3 

kg/m  which  corresponds  to  sea  level  pressure  and  a  temperature 
of  0°C.  Elliptic  loading  (C  =  4/ir)  was  assumed. 

5.1.2  Measurements 

The  conventional  method  of  determining  the  initial  vortex 
strength  has  been  to  average  the  measured  vortex  strengths  for 
early  vortex  detections.  Because  of  aircraft  noise  and  the  re¬ 
quirement  that  the  vortex  drift  past  the  antenna  to  complete  a 
measurement,  vortex  detections  before  10  seconds  produce  unreliable 
strengths.  Consequently,  the  average  strength  of  vortices  detected 
between  10  and  15  seconds  age  was  selected  as  a  reasonable  indica¬ 
tion  of  initial  vortex  strength.  An  additional  requirement  of 
reasonably  high  vortex  transport  velocity  (>2  m/sec)  was  imposed 
to  reduce  the  influence  of  vortex  decay  during  the  measurement. 

For  example,  the  measurement  of  a  20-meter  average  strength  lasts 
for  20  seconds  for  a  transport  velocity  of  2  m/s ec.  Table  9  shows 
the  results  of  averaging  the  vortex  strength  over  vortex  detections 
between  10  and  15  seconds.  The  average  strength,  standard  devia¬ 
tion,  and  ratio  of  standard  deviation  to  average  strength  are 
listed  for  ten  aircraft  types  and  three  averaging  radii.  One 
should  note  that  most  of  the  vortex  detections  (about  75  percent) 
used  are  first  vortices  because  they  are  the  first  to  reach  a 
MAVSS  antenna.  Some  difference  was  noted  between  the  first  and 
second  vortex  strengths,  the  latter  being  significantly  greater 
for  some  aircraft  types,  especially  for  5-m  averaging  radius. 

The  average  strengths  r’(10  m)  and  C’(20  m)  in  Table  9  were  used 
to  calculate  the  values  for  F'  (10-20  m)  in  Table  8  by  means  of 
Equation  5. 

Since  the  data  in  Section  2.1  show  indications  that  the  vortex 
decay  can  be  significant  even  in  15  seconds,  plots  were  generated 
showing  the  average  strength  of  vortices  detected  at  different 
ages.  It  was  hoped  that  the  strength  values  could  be  extrapolated 
back  to  zero  time.  Figures  11  through  13  are  samples  of  these 
plots;  the  plots  for  the  other  aircraft  types  are  shown  in 
Appendix  D.  The  2nd,  4th,  and  6th  points  in  these  plots  are  not 
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TABLE  9.  AVERAGE  INITIAL  VORTEX  STRENGTH  (mVsec):  All  Landing  Vortex  Detections  Between 
10  and  15  Seconds  With  Transport  Speeds  Above  2  Meters/Second. 


DECAY  OF  AVERAGE  VORTEX  STRENGTH  FOR  ALL  B-747  DETECTED  VORTICES:  AVERAGING 


FIGURE  12.  DECAY  OF  AVERAGE  VORTEX  STRENGTH  FOR  ALL  B-707  DETECTED  VORTICES:  AVERAGING 
RADIUS  =  S  m  (lower  line),  10  m  (middle  line),  and  20  m  (upper  line). 


statistically  independent,  but  are  simply  averages  of  the  adjacent 
points.  The  results  of  these  plots  were  surprising.  It  was  found 
that  the  vortex  strength  decay  could  be  approximated  by  a  linear 
decay  to  zero  at  a  particular  time  which  turned  out  to  be  120 
seconds  for  aircraft  with  four  wing-mounted  engines  and  the  B-727, 
and  140  seconds  for  other  aircraft.  For  some  aircraft  types 
(see  Figures  11  and  12)  the  linear  decay  curve  fit  reasonably  well 
even  back  to  zero  age.  For  others  (see  Figure  13)  the  vortex 
initially  decays  more  rapidly  until  it  reaches  the  linear  decay 
region.  The  vortex  decay  implications  of  these  plots  will  be 
deferred  to  Section  5.2.  The  analysis  here  makes  use  only  of  the 
zero  time  intercepts  of  the  linear  decay  curves.  Table  10  lists 
these  strength  values  and  compares  them  with  the  average  strength 
values  of  Table  9. 

The  vortex  strength  versus  averaging  radius  data  can  be  used 
to  derive  the  core  radius  according  to  the  procedure  described  in 
Section  2.1.  Figure  3  (Equation  6)  is  used  to  obtain  the  core 
radius  from  the  ratio  of  the  strengths  at  two  different  averaging 
radii.  The  results  of  this  analysis  for  three  radius  pairs  are 
listed  in  Table  10.  Using  Equation  3,  the  vortex  circulation  rB 
is  calculated  from  the  core  radius  rc  and  the  strength  P  at  the 
larger  radius.  If  the  vortex  model  fits  the  data  well,  the  core 
radius  and  circulation  will  be  the  same  for  all  three  radius 
pairs.  If  not,  the  best  selection  for  evaluating  is  the  10/20- 
meter  pair  since  it  is  least  affected  by  what  happens  at  the 
vortex  core.  The  10/20  values  of  T  are  listed  in  Table  8.  The 
calculated  values  of  core  radius  range  from  3.1  to  6.3  meters. 

The  wide-body  aircraft  values,  which  range  from  4.6  to  6.3  meters 
are  significantly  larger  than  the  other  aircraft  whose  values  lie 
between  3.1  and  4.8  meters.  Apart  from  this  well-defined  differ¬ 
ence,  there  seems  to  be  no  consistent  variation  of  core  radius 
with  aircraft  size.  One  would  suspect  that  the  MAVSS  itself  is 
limiting  the  core  size  to  perhaps  a  minimum  of  3  meters  because  of 
the  effects  discussed  in  Section  2.2.  Conventional  estimates  of 
core  size  are  significantly  smaller  than  those  in  Table  10, 
especially  for  the  small  aircraft  with  tail-mounted  engines. 
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TABLE  10.  ANALYSIS  OF  INITIAL  VORTEX  STRENGTHS 


1 

1 

INPUT  DATA 
Averaging 

Radius  . 

(b)  r'(l)  r<2) 

Radius 

Pair 

RADIUS-PAIR  ANALYSIS 

r  (1)  r  (2)  r»(l> 

(m)  (m)  (o2/sec) 

r.tt) 

(m^/sec) 

Aircraft 

Tyoe 

DC-9 

5 

51 

51 

5/10 

470 

371 

178 

- 139 — 

10 

93 

84 

10/20 

4.8 

4.0 

200 

161 

20 

136 

117 

5/20 

4.2 

3.4 

191 

155 

B-737 

5 

47 

52 

5/10 

4.6 

4.4 

198 

204 

10 

94 

100 

10/20 

4.8 

3.5 

203 

179 

20 

138 

135 

5/20 

4.7 

4.0 

201 

186 

8-727 

5 

60 

67 

5/20 

4.3 

3.8 

230 

222 

10 

115 

120 

10/20 

4.8 

4.0 

249 

230 

20 

169 

167 

5/20 

4.4 

3.8 

241 

228 

B-707 

5 

79 

86 

5/10 

3.5 

3.4 

239 

251 

l 

10 

136 

145 

10/20 

4.0 

4.0 

259 

278 

i 

20 

188 

202 

5/20 

3.6 

3.5 

251 

270 

1 

DC-8 

5 

85 

82 

5/10 

3.8 

3.8 

279 

270 

1 

10 

151 

146 

10/20 

4.2 

4.2 

298 

290 

j 

20 

213 

207 

5/20 

3.9 

3.9 

291 

294 

1 

DC-8H 

5 

84 

81 

5/10 

4.0 

4.0 

290 

281 

10 

152 

147 

10/20 

4.0 

4.0 

289 

281 

20 

210 

204 

5/20 

3.9 

3.9 

287 

286 

B-707H 

5 

83 

94 

5/10 

3.8 

3,5 

274 

283 

10 

148 

161 

10/20 

3.2 

3.5 

249 

286 

20 

193 

216 

5/20 

3.5 

3,4 

255 

286 

DC-10 

5 

61 

76 

5/10 

6.0 

4.4 

361 

301 

10 

138 

148 

10/20 

5.8 

5.2 

348 

341 

20 

218 

224 

5/20 

5.8 

4.7 

348 

332 

i! 

L-1011 

5 

69 

77 

5/10 

4.6 

4.9 

309 

359 

10 

147 

163 

10/20 

5.5 

5.2 

354 

377 

20 

227 

248 

5/20 

5.4 

5.2 

351 

374 

B-747 

5 

87 

102 

5/10 

5.0 

5.0 

405 

473 

10 

181 

211 

10/20 

6.3 

5.4 

497 

499 

20 

299 

323 

5/20 

5.7 

5.0 

473 

487 

(1)  Average  of  Detections  Between  10  and  15  Seconds. 

(2)  Linear  Extrapolation  to  Zero  Age. 
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Definitive  laser  Doppler  velocimeter  (LDV)  data  on  the  landing 
configuration  velocity  profile  is  available  (Ref.  5)  on  the  B-747. 
The  maximum  tangential  velocity  occurs  at  2.5-meter  radius. 

Because  the  circulation  profile  is  significantly  different  from 
that  (Equation  2)  used  in  the  core  radius  analysis  here,  the  peak 
velocity  core  radius  cannot  be  directly  compared  to  the  present 
results.  If  one  applies  the  analysis  of  Table  10  to  the  LDV 
vortex  strengths,  one  obtains  a  core  radius  of  3  to  3.5  meters 
at  age  25  seconds.  The  core  radius  appears  to  decrease  about 
2.5  meters  at  age  60  seconds.  Table  1  illustrates  the  factors 
which  may  be  needed  to  correct  the  average  circulation  values  in 
Table  9  for  the  error  in  measuring  rc<  The  value  of  F'(5  m)  is 
likely  to  be  a  factor  of  2  too  low  for  the  B-747.  This  is  a  large 
error  which  will  significantly  affect  the  hazard  duration  estimate 
for  small  following  aircraft.  The  errors  for  the  other  aircraft 
types  cannot  be  estimated  with  any  great  confidence.  Consequently, 
the  sensitivity  of  the  results  to  such  errors  must  be  considered 
in  subsequent  analyses.  Any  analysis  of  vortex  decay  which  com¬ 
pares  the  strength  of  aged  vortices  to  the  initial  measured  strength 
is  likely  to  be  unaffected  unless  the  core  size  changes  signifi¬ 
cantly  as  the  vortex  decays.  In  any  case,  the  value  of  is  not 
affected  by  the  core  radius  error. 

Figure  14  shows  log-log  plots  comparing  the  measured  (10-15 
second  age)  values  of  F^  to  the  two  calculated  values  in  Table  8 
which  represent  variable  and  fixed  airspeed,  respectively.  The 
fixed  airspeed  values  come  somewhat  closer  to  being  proportional 
to  the  measured  values,  although  the  scatter  is  considerable. 

The  measured  values  are  generally  lower  than  the  calculated  ones, 
as  one  might  expect  for  a  number  of  reasons,  such  as  vortex  decay 
and  landing  weights  lower  than  the  maximum  value  used  to  calculate 
r^.  The  effects  of  airspeed,  air  density,  and  wing  loading  in 
Equation  15  are  more  difficult  to  assess.  The  line  corresponding 
to  constant  wing  loading  is  shown  in  Figure  14.  The  effect  of 
landing  flap  deployment  on  is  likely  to  be  relatively  small 
since  the  flaps  are  roughly  equivalent  to  a  wing  of  constant  load¬ 
ing  but  of  shorter  span.  The  net  effect  is  probably  a  slight  in¬ 
crease  in  r  . 

00 
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FIGURE  14.  COMPARISON  OF  MEASURED  TO  CALCULATED  VORTEX  CIRCULATIONS 


Another  way  of  correlating  the  measured  circulation  values 
is  shown  in  Figure  15  which  shows  log-log  plots  of  .circulation 
versus  aircraft  wingspan.  The  10-15  second  age  circulations 
(Figure  15a)  appear  to  correlate  somewhat  better  with  wingspan 
than  with  the  calculated  values  in  Figure  14.  The  straight  line 
in  Figure  15a  corresponds  to  F^  =  7.2b  (m^/sec) .  Of  particular 
interest  for  subsequent  discussions  is  the  abnormally  low  10-15 
second  age  value  for  the  B-707H  in  all  the  comparisons.  The 
B-707H  had  the  fewest  cases  by  a  factor  of  two  of  all  the  air¬ 
craft  types  (see  Table  2). 

The  linear  decay  intercept  circulation  appears  to  have  a  non 
linear  relationship  with  wingspan  in  Figure  15b  where  a  line  with 
a  slope  of  1.25  gives  a  better  representation  of  the  data  than 
would  be  possible  with  a  slope  of  1.00.  The  same  relationship 
appears  to  hold  in  Figure  16  where  the  linear  decay  intercept 
circulation  is  compared  with  the  10-15  second  age  circulation. 

5.2  VORTEX  DECAY 

An  understanding  of  wake  vortex  decay  is  needed  to  evaluate 
wake  vortex  separation  criteria.  The  goal  of  the  present  effort 
has  been  to  find  a  simple  way  of  relating  the  decay  of  the  wake 
vortex  hazard  to  a  small  number  of  parameters  of  the  generating 
and  following  aircraft.  Sections  3.2  and  3.3  and  Appendix  B 
developed  methods  and  models  for  describing  the  decay  of  vortices 
from  a  particular  aircraft.  This  section  will  examine  how  the 
decay  depends  upon  aircraft  type. 

5.2.1  Average  Strength 

Figures  11-13  and  Appendix  D  show  that  the  decay  of  the 
average  strength  of  detectable  vortices  can  be  reasonably  fitted 
by  straight  lines  terminating  with  zero  strength  at  times  of  120 
or  140  seconds.  Although  the  selection  of  the  exact  termination 
time  and  initial  strength  is  somewhat  arbitrary,  nevertheless  the 
quality  of  the  straight  line  fits  is  remarkable.  A  mathematical 
interpretation  of  these  straight  lines  is  not  easy  to  define 
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(LINEAR  DECAY  INTERCEPT)  Tm  (10-15  SECOND  AGE)  (m  Vsec) 


WING  SPAN  (m) 


FIGURE  15.  DEPENDENCE  OF  MEASURED  VORTEX  CIRCULATION  UPON 
AIRCRAFT  WINGSPAN:  (a)  10-15  SECOND  AGE,  (b)  LINEAR  DECAY 
INTERCEPT 


(LINEAR  DECAY  INTERCEPT)  (m  /sec) 


FIGURE  16.  COMPARISON  OF  TWO  CIRCULATION  VALUES:  LINEAR  DECAY 
INTERCEPT  VERSUS  10-15  SECOND  AGE  AVERAGE 
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since  they  represent  only  those  vortices  with  strength  above  the 
MAVSS  detection  threshold,  whose  number  is  steadily  diminishing. 

The  following  "picture"  of  vortex  decay  is  a  plausable  explanation 
of  how  the  straight  lines  might  arise. 

The  straight  line  decay  for  all  averaging  radii  corresponds 
to  a  vortex  whose  velocity  profile  is  constant  but  whose  strength 
is  decreasing.  Table  10  evaluates  the  core  radii  for  the  straight 
line  decay.  In  some  sense  the  straight-line-decay  shape  represents 
a  vortex  "equilibrium"  decay  state.  For  some  cases  (Figures  11 
and  12),  the  "equilibrium"  decay  state  is  equivalent  to  the  initial 
state.  For  others  (Figure  13),  the  vortex  undergoes  a  rapid 
transition  from  the  initial  state  to  the  "equilibrium"  state. 

Table  10  shows  that,  in  general,  the  core  radius  is  somewhat 
smaller  for  the  "equilibrium"  state  than  for  the  initial  state. 

This  change  is  consistent  with  the  observation  (Ref.  5)  that  the 
vorticity  in  the  outer  portion  of  the  vortex  tends  to  diffuse 
away  more  rapidly  than  that  in  the  core.  It  is  appealing  to 
identify  the  linear  vortex  decay  with  turbulent  diffusion  which 
slowly  eats  away  at  a  vortex.  The  loss  of  vortices  below  the 
detection  threshold  would  then  be  identified  as  the  result  of 
instabilities  which  cause  a  vortex  to  rapidly  lose  strength  and 
become  undetectable.  The  inverse-time-squared  decay  in  the  models 
of  Appendix  B  corresponds  to  this  rapid  decay.  (Figures  7  and 
8  of  Volume  II  show  how  the  detection  rate  of  vortices  decays  in 
time)  . 

5.2.2  Hazard  Probability 

In  Section  5.1  of  Volume  II  of  this  report  it  was  found  that 
the  vortex  hazard  probabilities  for  different  aircraft  types  could 
be  characterized  by  assigning  an  effective  vortex  strength  r’eff 
to  a  generating  aircraft  such  that  the  vortex  hazard  probability 
can  be  expressed  as 

F(Ff,t)  =  Y (t  »rj/r^ff) ,  (16) 


5-17 


where  Y  is  a  universal  decay  function  for  all  aircraft  types.  One 
would  hope  that  the  variation  of  fgff  with  averaging  radius  can  be 
represented  by  an  effective  core  radius  rg££  which  might  be  the 
same  for  all  aircraft  types.  In  this  section  an  empirical  quest 
for  the  universal  function  Y  will  be  described.  The  next  section 
will  adopt  a  more  analytical  approach. 

The  first  choice  for  r^££  was  simply  the  initial  vortex 

strengths  (10-15  second  age)  F^  given  in  Table  9.  Figure  17  shows 

the  data  of  Figure  10  plotted  in  a  form  where  the  hazard  threshold 

r '  is  normalized  to  r'.  The  points  on  a  curve  in  Figure  17  cor- 
i  o  ? 

respond  to  the  value  of  (vortex  age/100)  where  the  curves  in 
Figure  10  pass  through  a  particular  value  of  hazard  probability  F. 
A  longer  time  is  required  to  reach  a  lower  value  of  F.  Likewise, 
the  time  to  reach  a  given  hazard  probability  is  longer  if  the 
hazard  strength  threshold  is  reduced. 

Figure  18  shows  how  the  averaging  radius  affects  the  plots 
of  Figure  17.  In  general,  the  three  averaging  radii  agree 
reasonably  well  in  such  normalized  plots;  changes  in  F^  of  20  per¬ 
cent  or  less  would  bring  the  curves  into  good  agreement.  The  data 
in  Figures  17  and  18  are  subject  to  distortions  below  of 

about  0.4  because  of  the  MAVSS  detection  threshold.  The  sharp 
break  in  the  curves  near  r^,/ r ^  =  1  occurs  when  becomes  higher 
than  some  of  the  initial  vortex  strengths.  The  portion  of  the 
curves  between  these  two  limits  can  be  represented  as  a  straight 
line  on  these  log-log  plots.  A  slope  of  -1.00  is  equivalent  to 
the  inverse  time  squared  vortex  decay  which  was  used  to  generate 
the  similar  theoretical  plots  in  Figure  37.  The  power  law  of  the 
vortex  decay  is  poorly  defined  in  these  plots  because  of  the  short 
range  of  the  data  validity. 

Figure  19  compares  the  normalized  decay  curves  for  all  the 
aircraft.  The  selected  hazard  probability  F  =  0.05  is  low  enough 
to  give  reasonable  trends  for  safety  considerations  but  is  still 
high  enough  to  have  reasonable  statistical  accuracy  for  most 
aircraft  types.  (The  B-707H  is  lost,  however,  because  of  too 
few  cases).  The  left  plots  of  Figure  19  show  the  curves  where 
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FIGURE  18.  TIME  TO  REACH  A  GIVEN  HAZARD  PROBABILITY  VERSUS 
NORMALIZED  HAZARD  THRESHOLD  T-f  / TA  FOR  THREE  AVERAGING  RADII: 

5  METERS  (NARROW  LINES),  10  METERS  (MEDIUM  LINES)  AND  20  METERS 
(BROAD  LINES) 


FIGURE  19.  COMPARISONS  OF  NORMALIZED  VORTEX  HAZARD  DECAY  CURVES  FOR  ALL  AIRCRAFT  TYPES 
F  =  0.05,  VORTEX  2,  AVERAGING  RADIUS  =  5,  10,  AND  20  METERS:  (a)  ri  =  INITIAL  VALUES, 
(b)  r'  ADJUSTED  TO  GIVE  GOOD  AGREEMENT  FOR  10  METER  AVERAGING  RADIUS 


was  taken  as  the  initial  10-15  second  mean  values.  These  values 
were  adjusted  to  give  excellent  agreement  in  the  middle  of  the 
-1.00  slope  region  for  10-meter  averaging.  The  required  changes 
in  were  -9  percent  for  DC-9,  B-727,  B-707,  and  DC-8H;  +  9 
percent  for  B-737;  +  20  percent  for  DC-10;  and  +  31  percent  for 
L-1011.  The  same  percentage  changes  were  made  for  all  radii. 

The  resulting  plots  overlay  fairly  well  for  all  radii.  The  only 
major  deviations  are  the  5-i.ieter  DC-8  and  DC-8H  curves.  Figure  20 
plots  these  selected  values  of  vortex  strength  as  an  effective 
reff  versus  aircraft  wingspan.  Lines  of  slope  1.00  are  drawn 
through  the  B-727  points  as  a  guide  for  comparison. 

This  empirical  method  of  searching  for  an  effective  vortex 
strength  is  not  as  fruitful  as  the  analytical  approach  to  be  dis¬ 
cussed  in  the  next  section.  The  concept  of  an  effective  strength 
to  describe  the  vortex  hazard  is  not  completely  satisfactory  since 
the  agreement  of  the  curves  in  Figure  19  could  probably  be  improved 
with  displacements  in  the  vertical  time  axis  rather  than  the 
horizontal  strength  axis.  In  the  vortex  decay  modeling  these  two 
are  related  by  Equation  35  for  a  given  hazard  probability. 

Changing  the  time  axis  rather  than  the  strength  axis  allows  a 
more  consistent  treatment  of  the  effects  of  the  distribution  of 
initial  vortex  strengths  (see  Section  B.3). 

5.2.3  Decay  Modeling 

Analytical  models  of  vortex  decay  are  more  useful  than  the 
empirical  fits  of  the  last  section  because  they  can  be  extrapolated 
beyond  the  limits  of  the  measurements  to  lower  values  of  vortex 
strength  and  lower  probabilities.  Such  an  extrapolation  is  needed 
to  reach  the  actual  hazard  probabilities  encountered  operationally 
and  to  explore  the  possible  range  of  the  fractional  roll  control 
parameter  f. 

5. 2. 3.1  Unsatisfactory  Models  -  Before  discussing  the  vortex 
models  which  can  be  used  to  represent  vortex  decay  in  a  practical 
way,  two  unrealistic  vortex  models  which  are  sometimes  used  to 
normalize  vortex  decay  are  examined.  The  first  is  the  similarity 
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WING  SPAN  (m) 


FIGURE  20.  EFFECTIVE  VORTEX  STRENGTH  FOR  DETERMINING  THE  DECAY 
OF  THE  WAKE  VORTEX  HAZARD 
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model  which  is  consistent  with  the  eddy  viscosity  decay  of 

vortices.  This  model  assumes  that  the  tangential  velocity  profile 

v  of  the  vortex  depends  upon  radius  r  and  age  t  in  the  form 
2 

v(r  /t).  In  this  model  the  vortex  core  radius  r  slowly  grows 

1/2  ^ 
as  t  '  and  the  value  of  remains  constant.  In  fact,  real 

vortices  tend  to  have  constant  or  decreasing  core  radii  with 

decreasing  r^. 

The  second  model,  which  is  used  to  normalize  scale-model  data, 
assumes  that  the  decay  of  wake  vortices  can  be  normalized  to  the 
number  of  wingspans  downstream  multiplied  by  several  correction 
factors : 


r' (t)/r^  =  z(tvcL/bA),  (17) 

where  V  is  the  airspeed,  is  the  lift  coefficient,  b  is  the 
wingspan,  A  is  the  wing  aspect  ratio,  and  Z  is  a  universal  decay 
function.  In  Sections  B.2  and  B.3  the  vortex  decay  time  is  found 
to  have  a  much  smaller  dependence  upon  wingspan  than  given  by 
Equation  17.  The  unsuitability  of  Equation  17  for  describing 
vortex  decay  in  the  atmosphere  may  be  at  least  partly  due  to  the 
fact  that  the  experimental  data  are  related  to  the  statistics  of 
vortex  decay  under  varying  conditions  while  Equation  17  describes 
the  deterministic  decay  of  a  single  vortex. 

5. 2. 3. 2  Simple  Stochastic  Model  -  Appendix  B  develops  a  variety 
of  analytic  models  which  describe  the  decay  of  wake  vortices.  One 
of  the  models,  the  simple  stochastic  model  (Section  B.3)  was 
selected  to  represent  vortex  decay  in  subsequent  discussions 
(Section  6.0).  The  form  of  that  model  is  outlined  here.  The 
example  examined  (DC-8)  will  be  of  particular  interest  in  Section 
6.0. 

The  vortex  decay  model  used  in  this  report  serves  the  function 
of  representing  the  vortex  hazard  decay  data  (as  in  Figures  8  to 
10)  in  an  analytical  form.  This  analytical  form  is  used  to 
evaluate  the  vortex  hazard  duration  for  values  of  hazard  probability 
F  and  hazard  threshold  Tj  which  lie  outside  the  region  where  direct 
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measurement  is  possible.  Only  two  free  parameters,  a  and  F',  are 
used  to  represent  the  vortex  hazard  decay  for  a  particular  air¬ 
craft  type  and  averaging  radius.  Figures  21  and  22  (which  cor¬ 
respond  to  Figures  38  and  39)  illustrate  this  model  for  the  case 
of  DC-8  vortex  2  with  15-meter  averaging  radius,  which  will  be 
important  for  the  discussion  in  Section  6.1,  Only  vortex  2  is 
used  since  the  results  will  be  used  to  evaluate  the  wake  vortex 
hazard  for  aircraft  using  the  same  runway.  Figure  21  shows  model 
fits  to  the  vortex  hazard  decay  using  a  separate  value  of  a  for 
each  hazard  threshold  r|  and  taking  the  10-second  value  of  hazard 
probability  as  the  initial  value.  From  Figure  21  a  value  of 

a  =  12.2  seconds  was  selected  to  fit  the  decay  curves  (r4-  =  75  and 

2  1 

100  m  /sec)  where  the  measurements  are  valid.  This  single  value 

of  a  is  used  to  generate  the  model  curves  for  all  values  of  Tj  in 
Figure  22.  The  initial  strength  distribution  assumes  a  20  percent 
standard  deviation  in  the  initial  vortex  strength  (T*  =  182  m  /sec). 
This  method  of  fitting  the  data  gives  hazard  duration  times  that 
are  longer  than  measured  for  low  threshold  Fj  values.  The  fitted 
form  is  thus  conservative  in  estimating  the  vortex  hazard.  It  is 
also  likely  to  be  a  better  estimate  of  vortex  decay  than  the 
actual  data  since  the  MAVSS  cannot  detect  weak  vortices.  (Note: 
the  fitted  curves  also  overestimate  the  duration  of  the  vortex 
hazard  for  Tj  >  ,  a  region  of  little  practical  consequence.) 

The  simple  stochastic  model  thus  is  a  realistic  representation 
of  the  decay  of  the  wake  vortex  hazard  for  the  purposes  of  the 
analysis  of  Section  6.0.  For  aircraft  types  and  averaging  radii 
other  than  the  case  considered  here,  the  value  of  comes  from 
Table  9  and  the  value  of  a  from  Figure  40. 


5-27 


IT/J00)**2  (5EC**2) 


FIGURE  21.  FITTED  PROBABILITY  DECAY  CURVES  FOR  THE  DC- 8, 
15-METER  AVERAGING  RADIUS,  AND  VORTEX  2.  THE  SIX  CURVES 
REPRESENT  T'  VALUES  OF  200,  150,  100,  75  ,  50,  AND  30  mVsec 
RESPECTIVELY,  FROM  LEFT  TO  RIGHT.  THE  MEASURED  DATA  ARE 
PLOTTED  WITH  THICKER  LINES  THAN  THE  FITTED  CURVES.  THE 
NUMBERS  BELOW  THE  CURVES  ARE  THE  VALUES  OF  j  (sec)  FOR  EACH 
CURVE . 


6.  R  E  C AT  E  GOR I Z  AT ION 


In  this  section  a  preliminary  look  at  recategorizing  air¬ 
craft  based  on  vortex  behavior  is  presented.  Resources  did  not 
permit  a  comprehensive  analysis;  however,  the  method  of  analysis 
is  described  through  the  development  of  a  four-category  system. 
Recommendations  for  further  work  are  proffered. 

The  philosophy  to  be  followed  in  examining  the  safety  of 
aircraft  separation  standards  is  similar  to  that  used  in  previous 
vortex  safety  analysis  (Ref*  11).  The  present  separation  stand¬ 
ards  are  assumed  to  be  safe  for  those  pairs  of  generator/follower 
aircraft  which  occur  frequently  in  normal  airport  operations. 

The  vortex  hazard  probability  F  found  acceptable  in  current 
operations  will  be  assumed  to  be  acceptable  under  revised  vortex 
categorization.  The  hazardous  vortex  encountered  rate  R^ve  is 
given  by 


Rhve 


R  _  H  F 
op 


(18) 


where  RQp  is  the  rate  of  operations,  H  is  the  probability  of  en¬ 
countering  a  vortex,  and  F  is  the  probability  that  the  vortex  is 
hazardous.  The  acceptable  level  for  F  can  be  surprisingly  large 
for  several  reasons.  First,  the  probability  H  of  encountering  a 
vortex  is  very  small  because  of  the  normal  motion  of  vortices. 
Second,  a  large  fraction  of  the  aircraft  operations  occur  at 
separations  larger  than  the  minimum.  Since  F  decreases  rapidly 
with  separation  distance,  these  operations  contribute  little  to 
the  hazardous  encounter  rate.  In  fact,  there  is  probably  a  cor¬ 
relation  between  the  vortex  hazard  duration  and  aircraft  separa¬ 
tion  which  further  reduces  R,  .  The  operations  rate  R  tends 

nve  op 

to  be  low  under  conditions  such  as  early  morning  when  the  vortex 
lifetime  is  longest.  The  considerations  in  this  report  will  be 
limited  to  vortex  encounters  occurring  on  final  approach  inside 
the  middle  marker.  This  location,  where  the  MAVSS  landing  data 
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were  collected,  is  most  hazardous  for  vortex  encounters  for  two 
reasons:  (1)  The  normal  descent  of  wake  vortices  is  arrested  by 

the  ground  so  that  a  vortex  (vortex  2)  can  remain  stalled  in  the 
normal  approach  path.  (2)  The  low  aircraft  altitude  leaves  little 
room  for  recovery  from  a  vortex  encounter.  Reference  11  examined 
the  wake  vortex  hazard  both  here  and  at  higher  altitudes. 

6.1  CURRENT  EXPOSURE  TO  VORTEX  HAZARDS 

The  current  IFR  runway- threshold  wake-vortex  separation 
standards  are  listed  in  Table  11.  Aircraft  are  separated  by 
certificated  maximum  gross  take-off  weight  W  into  three  categories 
Small  (W  <  12,500  lbs),  Large  (12,500  <  W  <  300,000  lbs),  and 
Heavy  (W  >  300,000  lbs).  The  spacings  of  3,  4,  5,  and  6  nautical 
miles  translate  into  time  separations  of  80,  107,  133,  and  160 
seconds,  respectively,  assuming  a  nominal  approach  speed  of  135 
knots.  The  most  notable  feature  of  the  current  separation  stand¬ 
ards  is  the  extremely  large  range  of  aircraft  sizes  within  the 
Large  category. 

The  procedure  for  fitting  the  vortex  decay  data  described  in 
Section  5.2.3  is  used  to  evaluate  the  vortex  hazard  probability  F 
at  the  times  (80,  107,  133,  or  160  seconds)  corresponding  to  the 
minimum  separations  in  Table  11.  Table  12  shows  the  results  for 
f  =  1.0  (Equation  13  and  Table  4),  and  Table  13  the  results  for 
f  =  0.5.  Different  averaging  radii  equal  to  the  following  air¬ 
craft  semispan  are  used  as  appropriate  to  the  size  of  the  aircraft 
The  hazard  probability  results  in  Tables  12  and  13  show  some 
departures  from  the  expectation  that  larger  aircraft  should 
generate  more  hazardous  vortices.  The  B-737  and  the  L-1011  show 
vortex  hazards  similar,  respectively,  to  the  larger  B-727  and 
B-747.  It  is  tempting  to  ascribe  this  increased  hazard  to  the 
greater  persistence  (see  Section  B.2)  of  vortices  from  aircraft 
with  two  wing-mounted  engines  (e.g. ,  DC-10,  L-1011,  B-737). 

The  highest  probabilities  in  Tables  12  and  13  occur  for  IO¬ 
meter  semispan  Large  (e.g.,  Gulfstream  II)  aircraft  behind  the 
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TABLE  11.  VORTEX  IFR  SEPARATION  STANDARDS  AT  RUNWAY  THRESHOLD 

(SEPARATIONS  IN  NAUTICAL  MILES) 


LEADING 

AIRCRAFT 

HEAVY 

LARGE 


HEAVY 


4 


FOLLOWING  AIRCRAFT 
I  LARGE 


SMALL 


5 


6 


3 


4 


SMALL 


3 


3 


3 


TABLE  12.  HAZARD  PROBABILITY  AT  CURRENT  MINIMUM  SEPARATION  FOR  f=1.0 


FOLLOWING  AIRCRAFT 

HEAVY 

LARGE 

SMALL 

SEMI SPAN 

20m 

20m 

15m 

10m 

10m 

5m 

LEADING 

AIRCRAFT 

HEAVY 

B-747 

2.2xl0~4 

L-1011 

9.9xl0"6 

1.8xl0~6 

DC- 10 

B-707H 

DC-8H 

LARGE 

DC-8 

2.4xl0~4 

2.4xl0*4 

1 . 3xl0~3 

1 . OxlO2 

2.0xl0_5 

2 . 3xl0~4 

B-707 

2.3xlO~5 

2.3xl0_5 

1.6xl0~4 

5 . Ixl0~3 

3.3xl0-6 

7.8xl0-5 

B-727 

4.2xl0~6 

4.8x10 

B-737 

1.2xl0-6 

2.9xl0_4 

8.9xl0_8 

DC- 9 

7 . 8xl0_8 

7.2xl0"5 

-II 


i 


) 

l 


Those  combinations  with  no  entry  have  a  probability  less  than  10 


TABLE  13.  HAZARD  PROBABILITY  AT  CURRENT  MINIMUM  SEPARATION  FOR  f=0.5 


FOLLOWING  AIRCRAFT 

HEAVY 

LARGE 

SMALL 

SEMISPAN 

20m 

20m 

15m 

10m 

10m 

5m 

LEADING 

AIRCRAFT 

HEAVY 

B-747 

2.3xL0'2 

6.9xl0~4 

1.4xl0~3 

7.1xlO~3 

1 . 8xl0~4 

6.4x10  4 

L-1011 

1.9xl0"2 

4.8x10  4 

1 . 3xlO-3 

l.lxlO-2 

3.9xlO~4 

2.3xlO~5 

DC-10 

5. 6x10  3 

5. 0xl0~5 

2.0xl0~4 

4.3xl0_3 

7.9xlO~5 

9.2xl0“5 

B-707H 

2. 3xl0_3 

9. Ixl0~6 

2 . 0xl0~4 

2.8xl0_3 

3.8xlO~5 

1.6x10  4 

DC-8H 

3.2xl0'4 

4.0xl0~5 

LARGE 

DC-8 

6.7xl0"2 

6. 7xl0~2 

l.lxlO-1 

2 . 3xlO_1 

2.2xlO~2 

2.  lxlO-2 

B-707 

3.3xl0"2 

3.3xl0~2 

5. 7xlO~2 

1.8xl0_1 

1 . 3xlO~2 

1.2xlO_2 

B-727 

7 . lxlO-3 

7 . Ixl0~3 

1 . 8xl0~2 

8 . lxlO-2 

2 . OxlO~3 

9.6x10  4 

B-737 

7. 3xl0_3 

7. 3x1 0-3 

2 . OxlO~2 

8 . 6xlO~2 

2 . 3xlO~3 

1.8xlO-3 

DC-9 

l.lxlO-3 

l.lxlO-3 

9. 7xlO~3 

5. 3xlO~2 

7 . 7xlO~4 

4. 2xl0~4 

Those  combinations  with  no  entry  have  a  probability  less  than  10  ^ 

*f=0.6  for  5-ra  averaging  radius. 

! 

i 
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Large  DC-8  and  B-707  aircraft.  This  combination  occurs  too  in¬ 
frequently  in  normal  operations  to  constitute  a  satisfactory 
choice  for  a  safe  vertex  hazard  probability.  The  next  highest 
probability  in  the  tables  occurs  for  a  15-meter  semispan  aircraft 
following  a  Large  DC-?  Since  this  size  follower  represents  the 
DC-9  and  the  B-737  ,  this  combination  occurs  often  enough  in  normal 
operations  to  represent  a  safe  level  of  hazard  probability.  With 
the  increase  in  commuter  traffic  since  deregulation,  the  number 
of  aircraft  operations  in  the  10-meter  semispan  Large  category 
has  increased  considerably.  This  higher  hazard  probability  for 
such  aircraft  following  the  biggest  Large  aircraft  raises  possible 
questions  concerning  the  safety  of  the  current  separation  stand¬ 
ards  . 

The  DC- 8  with  a  15-meter  semispan  following  aircraft  repre¬ 
sents  the  highest  vortex  hazard  probability  which  is  known  to 
be  safe.  Figure  23  shows  how  the  hazard  probability  for  this  case 
depends  upon  the  selected  value  of  f  ,  which  is  the  ratio  of  the 
maximum  safe  vortex- induced  rolling  moment  to  the  nominal  maximum 
control  -  induced  rolling  moment.  This  ratio  is  given  by  f  =  r,j,/5b 
according  to  Equations  13  and  14,  The  values  for  f  =  1.0  and  0.5 
were  included  in  Tables  12  and  13,  respectively. 

The  shape  of  the  F.  versus  f  curve  is  virtually  identical  to 

2  ^  5 

the  F  versus  t  plots  such  as  Figure  9.  This  similarity  stems 
from  the  fact  that  f  is  proportional  to  F2,  and  the  model  used 
gives  F  as  a  function  of  r,j,  t  as  long  as  r,j,  is  less  than  70  per¬ 
cent  of  T'  (which  is  true  for  all  the  hazard  cases  to  be  evaluated 

G 

in  the  following  discussion). 

The  relationship  between  F  and  fg  in  Figure  23  can  be  used 
to  investigate  the  sensitivity  of  the  safe  vortex  separation  time 
to  the  assumed  value  of  f  .  (Here  the  subscript  s  refers  to  the 
value  of  f  for  the  known  safe  case.)  For  example,  if  a  10-meter 
semi  span  aircraft  has  the  same  hazard  model  parameters  (p  and  K  in 
Equation  12)  as  the  DC-9  and  the  B-737,  then  the  hazard  duration 
can  be  calculated  as  the  time  for  the  10-meter  hazard  probability 
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to  reach  the  value  F  using  the  same  value  of  f  =  f  (r™  =  f5b) . 

The  resulting  safe  time  t  behind  a  DC-8  is  90.0  seconds  for  f 

s  s 

equal  to  0.7  or  less.  The  value  of  t  increases  to  90.8  seconds 

at  fg  =  1.0.  This  small  dependence  upon  f  is  a  consequence  of 

Tj  being  less  than  for  for  this  range  of  f  values.  The 

similar  functional  form  for  all  the  decay  curves  leads  to  very 

little  dependence  of  tg  upon  f  .  In  all  cases  to  be  presented 

here,  the  variation  in  t  for  fg  between  0.1  and  1.0  is  less  than 

2  seconds.  The  values  plotted  will  be  for  f  =  0.5.  The  value 

s 

of  t  =  90.0  seconds  for  a  10-meter  semispan  aircraft  following 
a  DC- 8  is  larger  than  the  defined  safe  time  of  80  seconds  for  a 
15-meter  semispan  aircraft  following  a  DC- 8,  as  one  would  expect. 
Tables  12  and  13  show  a  higher  hazard  probability  at  80  seconds 
for  the  smaller  semispan.  A  longer  time  must  elapse  before  the 
10-meter  hazard  probability  drops  to  the  80-second  value  for  15- 
meter  semispan. 

It  is  possible  that  the  values  of  p  and  K  for  the  10-meter 
semispan  aircraft  may  be  significantly  different  from  those  for 
the  15-meter  semispan  aircraft.  For  example,  the  value  of  p  for 
the  DC- 9  with  spoilers  is  known  (Ref.  12)  to  be  0.11  rather  than 
the  0.06  and  0.08  p  values  assigned  to  jet  transport  and  general 
aviation  aircraft,  respectively. 

Another  effect  which  affects  the  value  of  the  safe  time  t 

s 

in  the  same  way  as  differences  in  p  and  Y  is  the  error  introduced 
by  the  MAVSS  overestimate  of  the  core  radius.  According  to 
Table  1,  if  the  DC-8  core  radius  is  actually  1.5  m  instead  of  the 
measured  value  of  about  4.0  m,  the  10 -m  strength  will  read  a 
factor  of  1.31/1.50  =  0.87  low  relative  to  the  15-meter  refer¬ 
ence  value.  The  value  of  f  used  to  evaluate  the  safety  time  for 
a  different  semispan  should  be  reduced  by  this  same  factor.  A 
similar  reduction  in  f  is  needed  if  the  actual  value  of  p  is 
lower  than  the  0.11  DC- 9  value  with  spoilers;  a  reduction  factor 
of  0.08/. 1  1  =  .7  3  might  be  required  for  a  general  aviation  air¬ 
craft.  Because  of  the  possible  existence  of  these  correction 
factors,  an  additional  safe  time  ts  =  116.2  seconds  for 

6-8 


I 


* 


I 


f  =  0.6  f  is  calculated  to  -valuate  the  sensitivity  of  the 
results  to  these  two  types  of  effects.  The  time  is  significantly 
increased  since  these  effects  are  equivalent  to  reducing  the 
effective  hazard  threshold.  The  actual  hazard  duration  probably 
lies  between  the  f  =  fg  and  f  =  0.6  f  values  of  90.0  and  116.2 
seconds . 

Table  14  shows  a  similar  selection  of  variations  in  f/f 
which  are  needed  to  estimate  the  sensitivity  of  the  vortex  hazard 
duration  to  p  and  K  variations  and  to  core  size  errors  for  three 
aircraft  types  which  will  be  important  to  the  subsequent  discus¬ 
sion.  The  assumed  values  of  the  core  radii  are  the  values  con¬ 
tained  in  Table  1 . 

Since  the  DC-8  has  the  most  persistent  vortex  hazard  of  any 
of  the  Large  aircraft,  it  can  be  used  to  represent  the  safe 
separation  required  behind  the  Large  category.  The  dependence  of 
the  safe  separation  behind  Large  aircraft  upon  the  semispan  of 
the  following  aircraft  is  plotted  as  circles  in  Figure  24.  The 
solid  points  are  for  f/f  =  1  and  the  open  points  are  for  the 
values  of  f/f  given  in  Table  14.  The  B-747  similarly  has  the 
most  persistent  vortex  hazard  of  the  Heavy  aircraft  fexcept  for 
a  10-meter  semispan  follower  where  the  L-1G11  hazard  has  greater 
persistence).  The  B-747  safe  separations  (L-1011  for  10-meter 
semispan)  are  plotted  as  squares  in  Figure  24  to  represent  the 
Heavy  category.  Data  for  the  DC- 9  are  also  included  in  Figure  24 
to  illustrate  the  hazard  associated  with  the  middle  of  the  current 
Large  category.  The  selection  of  the  DC-8  and  B-747  to  represent 
the  worst  wake  generators  of  the  Large  and  Heavy  categories  could 
be  in  error  for  small  averaging  radii  if  the  f  core  correction 
of  another  aircraft  type  is  enough  larger  to  offset  the  apparent 
slower  decay  of  these  largest  aircraft  types. 

The  rei  alts  of  Figure  24  are  in  reasonable  agreement  with  the 
current  separation  standards  of  Table  11,  particularly  for  the 
Heavy  category.  The  5 -nautical -mi le  separation  for  Large  air¬ 
craft  behind  Heavies  is  satisfactory  for  all  semispans  for 
f/f  =  1  and  for  semi  spans  above  10.5  meters  when  considering  the 
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TABLE  14.  SENSITIVITY  ANALYSIS  VARIATIONS 


LEADING 

AIRCRAFT 

FOLLOWING 

AIRCRAFT 

SEMISPAN  (m) 

CORRECTION 

f/f5" 

DC-8 

15 

0 

1.00*** 

B-  747 

10 

0.87 

0.60 

DC- 9 

5 

0.60 

0.40 

B-  7  47 

20 

1 . 02 

0.70 

15 

0.95 

0.70 

10 

0.81 

0.60 

5 

0.  55 

0.40 

DC-9 

15 

0.90 

0.70 

10 

0.74 

0.60 

5 

0.44 

0.40 

*  Relative  to  15-meter  DC-8  case. 

ASSUMED  CORE  RADII: 

MEASURED  ACTUAL 

DC-8  4.0  1.5 

B-  747  5.5  2  .5 

DC-  9  4.0  0.5 

**  Incorporating  T '  core  corrections  and  0.73  factor  for  p  and  K. 
***  Defined  as  safe. 
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sensitivity  analysis  variations.  The  6-nautical -mi le  separation 
for  Small  aircraft  behind  Heavies  is  very  conservative  for  all 
semispans  for  f/f  =  1,  and  is  satisfactory  for  semispans  above 
6.5  meters  when  considering  the  sensitivity  analysis  variation. 
The  spacings  behind  the  Large  category  are  less  conservative  for 
the  assumed  variations  of  Table  14.  The  4-nautical-mile  separa¬ 
tion  for  Small  aircraft  is  satisfactory  for  all  semispans  for 
f/f  =  1  and  only  for  semispans  above  11  meters  when  considering 
the  variations  of  Table  14.  These  obser rations  set  the  stage  for 
a  proposed  four-category  separation  standard  to  be  examined  in 
the  next  section 

6.2  EVALUATION  OF  A  FOUR- CATEGORY  SYSTEM 

This  report  will  restrict  its  examination  of  possible  wake 
vortex  recategorization  to  a  single  scheme,  namely  dividing  the 
Large  category  into  two  categories  termed  Large  and  Medium.  The 
break  between  Large  and  Medium  is  set  just  below  the  DC- 9  and 
B-737  so  that  all  commercial  jet  transports  will  be  either  Heavy 
or  Large.  Such  a  split  is  in  use  in  the  United  Kingdom.  This 
split  can  simultaneously  improve  safety  and  reduce  traffic  delays 
because  it  more  accurately  reflects  the  actual  wake  vortex 
hazards.  A  glance  at  Table  11  shows  that  the  current  separations 
are  extremely  conservative  for  most  pairs  of  aircraft  and  that 
only  a  few  pairs  experience  any  significant  hazard.  The  use  of 
four  categories  rather  than  three  allows  a  more  uniform  hazard 
probabi 1 i ty . 

Table  15  shows  proposed  separation  standards  for  a  four- 
category  system.  Only  three  pairs  show  changes  from  current 
standards.  The  spacing  of  a  Large  behind  a  Heavy  is  reduced  from 
5  to  4  nautical  miles.  The  spacing  of  a  Medium  and  a  Small 
behind  a  Large  is  increased  from  3  to  4  and  from  4  to  5  nautical 
miles,  respectively.  The  net  effect  of  these  changes  is  likely 
to  be  an  increase  in  capacity  at  the  major  hub  airports  where  the 
traffic  is  predominatly  a  mix  of  the  Large  and  Heavy  aircraft  of 
the  four  -  category  system.  At  some  airports  the  Medium  and  Small 
aircraft  are  assigned  to  a  separate  shorter  runway. 
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TABLE  15.  PROPOSED  FOUR-CATEGORY  IFR  RUNWAY  THRESHOLD  SEPARATION 


(SEPARATIONS  IN  NAUTICAL  MILES) 


LEADING 

FOLLOWING  AIRC1 

1AFT 

AIRCRAFT 

HEAVY 

LARGE 

MEDIUM 

SMALL 

HEAVY 

4 

4*  * 

5 

6 

LARGE 

3 

3 

4* 

5* 

MEDIUM 

3 

3 

3 

4 

SMALL 

3 

3 

3 

3 

*  An  increase  from  current  separations. 
**  A  decrease  from  current  separations. 
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The  reduction  of  the  Large  behind  Heavy  spacing  from  5  to  4 
nautical  miles  is  well  justified  by  the  data  in  Figure  24.  The 
points  for  reduced  f/f  lie  only  slightly  above  the  4-nautical- 
mile  line  for  the  appropriate  15-  and  20-meter  semispans.  Most 
of  the  reduction  in  f/f  is  intended  to  account  for  p  Values  less 
than  those  of  the  DC- 9  and  B-737.  Since  those  are  exactly  the 
aixcraft  most  at  hazard  behind  a  Heavy  aircraft,  such  a  reduction 
is  unnecessarily  large;  the  actual  hazard  duration  is  therefore 
comfortably  below  the  4-nautical  mile  line. 

The  proper  safe  separation  of  Medium  and  Small  aircraft 
behind  Large  aircraft  is  less  easily  determined  because  of  the 
uncertainties  associated  with  the  MAVSS  core-size  errors.  A  strict 
requirement  that  the  separations  remain  below  the  open-circle 
line  in  Figure  24  leads  to  separations  which  seem  too  conservative 
in  comparison  to  the  current  standards.  The  4-nautical-mile 
Medium-behind-Large  separation  is  strictly  valid  only  for  follow¬ 
ing-aircraft  semispans  greater  than  11  meters.  Similarly,  the 
5 -naut ical -mi le  Small-behind-Large  separation  is  strictly  valid 
only  for  semispans  greater  than  7.5  meters.  It  is  apparent  that 
the  hazard  model  of  Section  3.1  related  the  vortex  hazard  for  a 
following  aircraft  to  the  wingspan  rather  than  the  weight  which 
is  used  to  define  the  current  categories.  For  the  purposes  of 
the  present  discussion  it  therefore  makes  sense  to  define  the 
boundaries  dividing  Large,  Medium,  and  Small  on  the  basis  of  wing¬ 
span.  The  Medium-Large  division  is  set  at  a  25-meter  wingspan 
(12.5-m  semispan)  which  is  just  below  the  DC-9  and  B-737.  The 
Small-Medium  division  is  more  arbitrary  and  is  set  at  a  15-meter 
wingspan  (7.5-m  semispan). 

The  proper  separation  of  a  Small  behind  a  Medium  is  difficult 
to  define  for  two  reasons.  First,  we  have  little  vortex  data  on 
Medium  aircraft.  Second,  the  aircraft  wingspan,  which  properly 
classifies  an  aircraft  as  a  follower,  may  not  properly  classify 
an  aircraft  as  a  generator.  The  effective  vortex  strength,  shown 
in  Figure  20,  is  roughly  proportional  to  wingspan  for  jet  trans¬ 
port  aircraft  which  have  similar  wing  loading.  However,  the 
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vortex  hazard  is  likely  to  be  much  less  for  an  aircraft  such  «s 
the  DC- 3  which  has  much  lower  wing  loading.  The  vortex  hazard 
from  the  DC-9,  which  is  just  above  the  Medium  category,  can  give 
some  clue  to  the  hazard  from  Medium  leading  aircraft.  The  DC-9 
data  in  Figure  24  show  that  4-nautical-mile  spacing  is  strictly 
safe  for  aircraft  with  semispans  above  7.5  meters.  In  accord 
with  the  philosophy  used  with  Large  leading  aircraft,  the  Small- 
behind-Medium  safe  separation  was  set  a  4  nautical  miles. 

6.3  HhAVY- LARGE  BOUNDARY 

The  data  in  this  report  have  some  bearing  on  current  ques¬ 
tions  regarding  the  division  between  the  Heavy  and  Large  cate¬ 
gories.  This  division  becomes  somewhat  less  important  in  the 
four-category  system  of  Section  6.2  than  in  the  current  system 
since  the  Large  and  Heavy  categories  have  the  same  separations  as 
followers . 

Volume  II  of  this  report  (Ref.  2)  examined  the  feasibility 
of  including  the  B-707H  and  DC-8H  in  the  Large  category.  The 
impact  of  this  change  on  the  Large  separations  in  Figure  24  is 
shown  in  Table  16.  The  DC-8H  actually  has  a  less  persistent 
vortex  hazard  than  the  DC- 8  and  is  not  included.  The  B-707H, 
however,  shows  a  more  persistent  hazard.  Because  of  the  small 
number  of  cases,  the  B-707H  results  have  relatively  poor  statis¬ 
tical  accuracy. 

A  number  of  new  jet-transport  aircraft  (e.g.,  A-300,  A-310, 
B-757,  B-767)  have  been  and  are  being  added  to  the  jet  transport 
fleet  since  the  MAVSS  data  were  collected.  They  all  have 
weights  near  the  break  point  between  the  Large  and  Heavy  cate¬ 
gories  and  they  all  have  two  wing-mounted  engines.  In  the  light 
of  the  observed  abnormally  high  persistence  of  vortices  from  air¬ 
craft  with  two  wing-mounted  engines,  it  may  be  advisable  to  assign 
such  aircraft  to  the  Heavy  category.  Such  has  been  the  practice 
in  the  UK  where  the  A-300  has  been  designated  a  Heavy  on  the  basis 
of  vortex  encounter  reports. 


TABLE  16.  COMPARISON  OF  SAFE  SEPARATIONS  FOR  THE 
DC-8  AND  THE  B-707H 


SEMI SPAN  OF 
FOLLOWING 
AIRCRAFT 
(m) 


DC-8 

B-707H 

t. 

t. 

f/fs 

s 

(sec) 

(sic) 

1.0 

80.0 

85.7 

1.0 

90.0 

97.0 

0.6 

116.2 

125.3 

1.0 

98.1 

112.0 

0.4 

155.0 

17  7.1 

6.4  RECOMMENDATIONS 

This  preliminary  analysis  indicates  that  a  four-category 
system  might  have  merit.  If  the  system  does  indeed  have  opera¬ 
tional  u+:lity,  a  host  of  questions  needs  to  be  examined  before 
proceding  <.0  implementation.  The  following  paragraphs  outline 
recommended  tasks. 

6.4.1  Safety  Analysis 

The  safety  analysis  herein  should  be  refined  by: 

1)  Investigating  the  adequacy  of  the  hazard  model  for  wide 
variations  in  aircraft  size.  (Should  a  pilot  response 
time  be  included?) 

2)  Examining  the  adequacy  for  the  MAVSS  measurements  (par¬ 
ticularly  the  core  corrections). 

3)  Combining  the  hazard  probabilities  for  each  aircraft  pair 
with  the  frequency  of  occurrence  for  each  pair  to  derive 
an  estimate  of  the  total  hazardous  encounter  rate  for 
the  current  and  proposed  systems. 

4)  Considering  encounters  at  locations  other  than  between 
the  middle  marker  and  touchdown  (including  takeoffs). 

6.4.2  Category  Boundaries 

The  selection  of  category  boundaries  (Heavy,  Large,  Medium, 
and  Small)  must  be  based  on  a  detailed  evaluation  of  the  exact 
characteristics  (e.g.,  p,  b,  and  W)  of  the  aircraft  which  lie  near 
the  boundary.  The  selection  of  boundaries  and  separations  must 
also  include  the  effect  on  capacity  at  different  types  of  airports 
such  as  commercial  hubs  and  busy  GA  fields. 

6.4.3  Operational  Questions 

The  air  traffic  implications  of  increasing  the  number  of 
categories  must  be  examined.  It  may  be  feasible  to  retain  a  three 
by  three  separation  matrix  but  with  different  leader  and  follower 
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categories  ( e.g .  ,  an  A-  300  ,  say,  be  classified  as  a  Heavy  when  it 
is  .he  leading  aircraft  but  classified  as  a  Large  when  it  is  the 
following  aircraft);  would  such  a  differentiation  be  useful? 

6.4.4  United  Kingdom  Experience 

The  UK  has  considerable  experience  with  both  operational  and 
safety  aspects  of  possible  recategorization.  Closer  cooperation 
with  the  UK  would  be  advantageous  to  answering  questions  and  also 
would  expedite  the  ICAO  acceptance  of  a  new  categorization  system. 
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APPENDIX  A 

MAVSS  DETECTION  THRESHOLD 


The  detection  of  wake  vortices  as  they  pass  over  successive 
MAVSS  antennas  is  accomplished  by  correlating  the  vertical  veloc¬ 
ity  profile  with  a  function  which  is  similar  to  a  vortex  signa¬ 
ture  (see  Volume  I,  Section  S.2).  For  vortices  from  landing 
aircraft  the  correlator  sums  vertical  velocities  over  a  region 
18  m  high  (5  range  gates)  and  20  m  wide.  The  threshold  assigned 
to  the  correlator  corresponds  to  a  uniform  absolute  velocity, 

Vc  =  2.5  m/sec  over  the  correlation  region.  Of  course,  the  verti¬ 
cal  velocity  signature  of  a  vortex  is  not  uniform  (see  Figure  1); 
it  is  zero  along  a  vertical  line  through  the  vortex  center,  reaches 
a  maximum  on  either  side  of  the  vortex  core,  and  then  decreases 
for  both  horizontal  and  vertical  displacements  from  the  vortex 
center  (Reference  13  shows  many  such  profiles).  Some  indication 
of  the  vertical  variation  in  average  velocity  can  be  obtained  by 
comparing  the  average  velocity  Vc  (absolute  value)  over  five 
range  gates  from  the  correlator  with  the  average  velocity  V^, 
along  the  single  range  gate  nearest  the  vortex  height.  The  ratio 
of  Vr  to  Vc  varied  between  l.1  .d  2.0  with  a  typical  value  of 
1.5  for  B-707  and  DC-8  data. 

One  estimate  of  the  vortex  strength  detection  threshold  can 
be  obtained  by  calculating  the  average  circulation  of  a  vortex 
with  a  constant  tangential  velocity  vQ.  According  to  equation  3, 
the  value  of  r'(r)  is  proportional  to  the  averaging  radius  r: 

r  '  (r)  =  TTVQr .  (19) 

If  we  take  vQ  =  2.5  m/sec,  the  average  velocity  Vc  of  the  corre¬ 
lator  threshold,  we  obtain  a  vortex  strength  threshold  estimate  of 

F’(r)  =  7.9  r.  (m2/sec)  (20) 
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For  consistency  this  value  should  be  reduced  by  the  ratio  Vr/Vc 
since  is  the  average  tangential  velocity.  Taking  Vr/Vc  =1.5 
leads  to  a  threshold  estimate  of 

r'  (r)  =  5.2  r.  (m2/sec)  (21) 

This  value  corresponds  to  the  hazard  fraction  value  f  =0.52  in 
Equation  13. 

A  more  direct  determination  of  the  effective  MAVSS  detection 

threshold  can  be  obtained  by  comparing  the  measured  strength  with 

the  value  of  the  detection  correlation.  Figure  25  shows  such  a 

comparison  for  10-m  averaging  radius  for  a  sample  of  B-707  and 

DC- 8  data.  To  increase  the  number  of  low  value  points,  data  for 

parts  of  the  sample  were  plotted  only  if  the  correlator  was  less 

2 

than  60  or  the  strength  less  than  80  m  /sec.  According  to  this 

2 

set  of  data,  vortices  of  strength  less  than  50  m  /sec  start  to  be 

undetected.  The  strength  for  a  50  percent  detection  probability 

2 

appears  to  be  about  30  m  /sec.  Note  that  the  10-m  averaging 
radius  was  used  in  Figure  25  because  the  correlator  and  strength 
are  averaged  over  the  same  distance  and  therefore  might  be 
expected  to  show  the  best  correlation. 

Another  way  of  estimating  the  detection  threshold  is  to 
examine  the  distribution  of  the  detected  vortex  strengths,  which 
is  shown  in  Figure  26  for  four  different  sizes  of  aircraft.  One 
would  expect  a  break  in  the  distribution  curve  below  the  detection 
threshold.  In  fact,  a  well-defined  break  is  observed  in  most 
cases.  The  dashed  lines  indicate  possible  continuations  of  the 
actual  vortex- st rength  distribution.  The  detection  thresholds 
appear  to  be  about  30  m  /sec  for  5-m  averaging  radius,  40  m  /sec 

7 

for  10-m  radius,  and  50  to  60  m “ / s e c  for  20 -m  radius.  These 

2 

numbers  were  used  (60  m^/sec  for  20-m  radius)  to  calculate  the 
limiting  values  of  f  in  Section  4.1. 
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FIGURE  25.  VORTEX  STRENGTH  VERSUS  AVERAGE  CORRELATION 
VELOCITY  Vc  FOR  A  SAMPLE  OF  B-707  AND  DC- 8  VORTICES 
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APPENDIX  B 
VORTEX  DECAY  MODELS 


An  accurate  model  for  vortex  decay  can  serve  a  number  of 
useful  purposes.  First,  it  would  allow  the  extrapolation  of 
decay  properties  beyond  the  experimental  limits  of  the  measure¬ 
ments.  Two  limits  are  of  particular  importance:  (1)  the  finite 
number  of  measurements  precludes  the  direct  measurement  of  very 
small  probabilities.  (2)  The  finite  MAVSS  detection  threshold 
distorts  the  vortex  decay  measurements  at  low  strength.  A  second 
purpose  for  a  decay  model  is  to  characterize  the  decay  by  a  small 
number  of  parameters  whose  variation  with  the  cases  selected  can 
be  used  to  determine  the  dependence  of  vortex  decay  upon  weather 
conditions  and  aircraft  type. 

B.l  SIMPLE  ANALYTICAL  MODEL 

Equation  8  in  Section  3.2  can  be  used  to  describe  the  long¬ 
time  asymptotic  decay  of  the  wake  vortex  hazard.  The  two  param¬ 
eters  of  the  model  a  and  d  represent,  respectively,  the  asymp¬ 
totic  slope  and  the  delay  of  the  decay.  The  delay  d  is  close  to 
zero  when  the  hazard  threshold  is  near  the  most  probable 
initial  strength  r^.  The  value  of  d  becomes  increasingly  positive 
as  rj.  decreases  below  F^  and  increasingly  negative  as  increases 

above  T'.  The  value  of  a  increases  as  rl  decreases  for  Ti  <  f 
o  1  1  o 

but  tends  to  remain  constant  for  Early  in  the  analysis 

of  the  data  a  simple  dependence  of  a  upon  F,j,  was  found,  namely 
that  the  slope  a  is  proportioned  to  F,j,  for  and  greater 

than  the  detection  threshold: 


r  =  c  r|.  (22) 

The  analysis  leading  to  this  result  for  the  data  in  Figure  5 

2 

is  presented  in  Table  17.  The  values  of  (t/100)  which  correspond 
to  a  probability  F(r.j,,t)  of  5  percent  will  be  used.  This  value  of 
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TABLE  17.  PARAMETERS  FOR  FIT  TO  DATA  IN  FIGURE  5 


r 

2  0 
fm  /sec) 

(T/100) 2 
for  F=0 . 05 

i\j.(T/ioo)2 
for  F=0 . 05 

30 

0.81 

24 

50 

0.68 

34 

75 

0.51 

38 

100 

0.  39 

39 

150 

0.09 

14 
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F  is  small  enough  that  the  delay  d  can  be  ignored  since  most  of 
the  variation  in  (t/100)  is  due  to  a  for  Fy  <  F^.  It  is  also 
large  enough  to  be  above  the  statistical  fluctuations  near  the 


single-case  level.  Use  Equation  8  with  d  =  0  to  relate  F  =  5 
percent  and  a  by  taking  the  natural  logarithm  of  both  sides: 

-  3.00  =  -  at2.  (23) 

Combining  Equations  22  and  23  yields 

G  =  3.00/(r|  t2).  (24) 

Table  17  shows  the  values  of  Fj(t/100)2  for  Figure  5.  These 


values  are  roughly  constant  in  the  region  between  the  detection 

2 

threshold  of  40  m  /sec  and  the  most  probable  initial  strength 

?  ~  fa 

f  of  135  rn  /sec.  The  corresponding  value  of  G  is  7.9  x  10* 

-2-1 
m  sec  . 

Although  the  assumption  of  d  =  0  in  the  analysis  of  Equations 

22  to  24  is  not  strictly  applicable  when  F-j,  is  significantly 

different  from  F^ ,  it  can  deal  with  many  of  the  cases  of  interest. 

Specifically,  the  values  of  T 'Q  in  Table  3  are  generally  close  to 

the  f  =  1.0  hazard  thresholds  of  Table  5.  The  ratio  of  F'  to  T™ 

o  I 

lies  in  the  range  of  0.8  to  1.85. 

The  assumption  d  =  0  is  Equation  8  combined  with  Equation  22 
leads  to  a  very  simple  analytical  expression  for  the  distribution 
function  P(r',t)  in  Equation  11: 

F(r^,t)  =  cxp[-GF|  t2]  =  /  P  (F '  ,t)  dF'  .  (25) 

rT 

Differentiation  with  respect  to  leads  to 

p(r’,t)  =  G  r2  cxp[-Gr|  t2].  (26) 
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The  decay  of  average  strength  from  Equation  10  becomes: 


r*  (t) 


(27) 


The  ensemble  mean  strength  decays  inversely  with  vortex  age 
squared.  One  should  note  that  this  model  describes  the  decay  of 
wake  vortices  with  a  single  parameter  G  which  depends  upon  air¬ 
craft  type,  meteorological  conditions,  and  strength  averaging 
radius . 

Figures  27  and  28,  respectively,  show  a  comparison  of  the 

predictions  of  Equations  25  and  26  with  the  measured  values  for 

all  B-707  vortices  and  a  10-m  averaging  radius.  In  Figure  27  the 

exponential  decay  lines  of  Equation  25  agree  exactly  at  F  =  5 

percent  for  =  75  and  100  m  /sec  because  those  points  were  used 

to  calculate  the  value  of  G.  The  straight-line  fit  to  the  100 

m  /sec  curve  (the  f  =  1.0  hazard  threshold)  agrees  very  well  for 

F  >  5  percent  but  agrees  less  well  for  lower  probabilities 

because  of  statistical  fluctuations.  The  agreement  is  poorer  for 

rj  =  75  and  50  because  the  assumption  of  d  =  0  is  less  reasonable. 

Figure  28  is  a  logarithmic  plot  of  Figure  6.  The  model  line 

2 

agrees  reasonably  well  at  70  seconds  between  r'  of  135  m  /sec 

2  ° 

and  the  detection  threshold  of  40  m  /sec,  as  anticipated  in  the 

evaluation  of  G.  For  earlier  times  (e.g.  ,  50  seconds)  where  the 

probability  distribution  has  a  well-defined  peak  the  fit  is  less 

good.  This  peak  in  the  probability  distribution  is  closely 

2 

related  to  the  need  for  nonzero  d  for  the  75  and  50  m  /sec  curves 
in  Figure  2’’ .  The  shape  of  the  model  curves  is  suitable  for  times 
of  90  and  110  seconds,  but  the  levels  are  low  because  of  statis¬ 
tical  fluctuations  and  the  loss  of  cases  below  the  detection 
threshold  of  40  m  /sec. 


B . 2  STOCHASTIC  MODEL 

The  simple  analytical  model  of  Section  B.l  suffers  from  a 
number  of  difficiencies ,  the  most  important  of  which  is  the  lack 
of  dealing  adequately  with  the  delay  in  the  vortex  decay.  A 


(VORTEX  AGE  (SEC)/ 100)  SQUARED 


FIGURE  27.  COMPARISONS  OF  MODEL  HAZARD  PROBABILITY  (STRAIGHT 
LINES)  WITH  MEASUREMENTS  (LINES  WITH  POINTS)  FOR  ALL  B^707 
VORTICES  WITH  10 -METER  AVERAGING.  THE  VALUES  OF  r»  (mZ /SEC) 
ARE  INDICATED.  T 


PROBABILITY  DISTRIBUTION  (SEC/m  )  X  10 


FIGURE  28.  COMPARISONS  OF  MODEL  PROBABILITY  DISTRIBUTIONS 
(STRAIGHT  LINES)  WITH  MEASUREMENTS  (LINES  WITH  POINTS)  FOR  ALL 
B- 707  VORTICES  WITH  10 -METER  AVERAGING.  THE  TIMES  (SEC)  AT 
WHICH  THE  DISTRIBUTION  IS  EVALUATED  ARE  INDICATED. 
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secondary  problem  is  the  lack  of  any  systematic  cons iderat i on  of 
the  initial  vortex  strength  distribution.  The  stochastic  model 
described  in  this  section  deals  directly  with  these  questions. 

The  basic  concept  of  the  stochastic  model  is  based  on  the 
observation  that  vortices  often  decay  in  two  stages.  During  the 
first  stage  the  vortex  strength  remains  roughly  constant.  during 
the  second  stage  the  strength  decays  rapidly.  This  type  of  decay 
can  be  modeled  mathematically  by 


r*  -  r  ( tl/t)n  t  >  tj 


r ' 


(28) 


where  r'  is  the  average  circulation  for  a  particular  aircraft 
and  averaging  radius,  F^  is  the  initial  value  of  F',  and  t^  is 
the  time  t  when  the  strong  vortex  decay  begins.  The  vortex  decay 
in  time  is  governed  by  the  power  n.  The  model  in  Section  B.l 
corresponds  to  n  =  2  (see  Equation  27) .  The  stochastic  feature 
of  the  model  is  introduced  by  assuming  that  the  time  the  decay 
starts  (tj)  is  distributed  in  a  normal  distribution  with  mean 
value  t  and  standard  deviation  a: 


P(t  )  =  J—  expT-ftj-t  )2  /2a2].  (29) 

/2tt  a  L  J 


The  three  parameters  of  the  model  are  thus  n,  t  ,  and  a. 

o 

If  the  distribution  in  initial  vortex  strengths  is  ignored, 
the  hazard  probability  F(r-j,,t)  can  be  calculated  explicitly  in 
terms  of  the  error  function.  The  limiting  value  of  initial 
decay  time  1 is  when  the  initial  strength  is  reduced  to 
according  to  Equation  28: 

t1T  =  t(F|/r^)1/n.  (30) 
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The  hazard  probability  is  given  by  the  integral  of  P(t^)  for  all 
times  greater  than  t ^ : 

00 

FCr^.t)  =  /  p Ctp  dtj  .  (31) 

"tlT 

This  integral  is  related  to  the  error  function  erf  which  satisfies 
the  relationship: 


7  r00  - 1 L 

1  -  erf (z)  «  —  J  e  dt .  (32) 

/T  z 

Combining  Equations  29  through  32  leads  to  the  expression 

F(I\J.,t)  =  \  [l  -  erf  ((t1T-t0)//2  a)l  .  (33) 

The  results  of  this  model  assuming  n  *  2  and  tQ  =  3o  are 
shown  in  Figure  29.  The  values  selected  for  r^/r^  correspond  to 
the  B-707  10-m  results  plotted  in  Figures  5  and  27.  The  curves 

7 

corresponding  to  =  50,  75,  and  100  m  /sec  are  in  reasonable 
correspondence  with  the  measurements  for  low  probabilities.  The 
30-m  /sec  curve,  however,  is  off,  perhaps  because  it  is  below  the 
MAVSS  detection  threshold  or  perhaps  because  of  the  wrong  value 
for  n. 

Several  features  of  Equation  33  and  Figure  29  are  worth 
noting.  First,  the  shapes  of  the  curves  are  similar  and  corre¬ 
spond  simply  to  scaling  the  horizontal  axis  for  different  values 
of  rj/T^.  One  way  of  looking  at  the  dependence  of  F  upon  Pj/Pq 
is  to  note  that  a  given  value  of  F  corresponds  to  a  particular 
value  of  the  argument  z  of  erf  in  Equation  33: 

z  =  [t(r^./r;)1/n  -  to]//2  o.  (34) 
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FIGURE  29.  VORTEX  DECAY  FOR  THE  STOCHASTIC  MODEL  WITH 

t  =  3o  AND  n  =  2 
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Equation  34  can  be  rearranged  to  give: 

1 

t/o  =  (t 0/o  +  n  z)  (r^/r^)  n  .  (35) 

A  second  feature  of  the  curves  in  Figure  29  is  that  the 
asymptotic  line  for  low  probabilities  is  not  straight  but  in  fact 
is  slightly  curved.  This  curvature  is  consistent  with  the 
asymptotic  expression  for  the  error  function: 

i  2 

(  1  -  erf(z))  -  — —  e'z  (36) 

/T  z 

The  inverse  z  dependence  in  Equation  36  causes  the  curvature. 

Equation  35  can  be  used  to  calculate  the  time  tp  required  to 
reach  a  certain  hazard  probability.  Listed  below  are  the  error- 
function  argument  z  values  leading  to  specified  hazard  probabil¬ 
ities,  F: 


F  z 


0.2 

0.59 

0.1 

0.91 

0.05 

1.16 

0.02 

1.45 

0.01 

1.65 

0.005 

1.82 

0.002 

2.03 

2 

Figure  30  shows  a  log- log  plot  of  the  dependence  of  (tp /a) 
upon  fj/fQ  for  three  F  values  for  the  choice  of  parameters 
tQ/c  =  3  and  n  =  2.  For  n  =  2  the  lines  have  a  slope  of  -1. 
Changing  the  value  of  n  simply  changes  the  slope  of  the  lines 
passing  through  the  point  where  =  1.  Changing  the  value  of 

tQ/o  moves  the  lines  up  or  down. 

The  stochastic  model  was  fitted  to  the  measured  data  on 
vortex  decay.  A  weighted  least-squares  fit  was  designed  to  find 
parameters  which  would  produce  reliable  predictions  of  the  vortex 
hazard  probabilities  at  long  times.  The  vortex  strength 
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distribution  is  evaluated  in  bins  of  width  20  m  /sec  (as  in 
Figure  6).  The  initial  strength  distribution  is  taken  as  the 
10-second  distribution.  The  stochastic  model  (parameters  a,  tQ, 
n)  is  used  to  calculate  the  probability  of  going  from  an  initial 
strength  bin  i  to  another  strength  bin  j  at  a  time  t  later.  The 
total  strength  distribution  at  time  t  is  obtained  by  summing  over 
the  initial  strength  distribution  P(i,10).  The  calculated 
strength  distribution  Pc(j,t)  was  compared  to  the  measured  dis¬ 
tribution  Pm(j,t)  by  evaluating  the  weighted  sum: 


t  N 

s=  £  E  (£n  P  (j ,t)  -  £n  P  (j,t))Z  N  (j)  (37) 

t=tx  j=l  c  m  m 

where  the  sum  is  taken  over  both  time  (even  multiples  of  10 

seconds)  and  strength  bins.  The  sum  S  is  minimized  to  yield  the 

least-square  fit.  The  use  of  logarithms  in  Equation  37  gives 

equal  weight  to  fractional  errors  in  the  distribution  at  both 

high  and  low  values.  The  sum  is  weighted  by  the  number  of 

measured  points  Nm(j)  in  a  bin  to  account  for  the  expected 

statistical  fluctuations.  Some  limitations  on  the  sum  are  placed 

to  achieve  the  goal  of  reliable  long-time  predictions.  First, 

the  lowest  strength  bin  is  set  to  be  above  the  MAVSS  detection 

2 

threshold.  The  lowest  bin  used  is  26  to  45  m  /sec  for  5-meter 

2 

averages,  46  to  65  m  /sec  for  10-meter  averages,  and  66  to  85 
2 

m  /sec  for  15-  and  20-meter  averages.  Second,  the  first  time 

used  is  40  seconds.  The  procedure  for  minimizing  S  is  to  select 

the  value  for  n  (1.0,  1.5,  2.0  or  2.5)  and  vary  a  and  t  until 

the  smallest  S  is  found.  The  smallest  increment  on  a  and  t  is 

o 

1/4  second. 

Figures  31  and  32  show  sample  results  of  the  least-square  fit 
for  two  aircraft  types,  B-727  and  the  B-747.  The  rms  error  in 
the  figures  is  the  root -mean- square  of  the  logarithm  difference 
in  Equation  37.  It  is  obtained  by  normalizing  S  to  the  sum  where 
the  quantity  squared  is  set  equal  to  1,  and  then  taking  the 


4 
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FIGURE  31.  LEAST-SQUARE  FIT  TO  VORTEX  DECAY :  ALL  B-727  VORTICES 


FIGURE  32.  LEAST-SQUARE  FIT  TO  VORTEX  DECAY:  ALL  B-747  VORTICES 


square  root.  An  rms  error  of  0.1  corresponds  to  an  rms  11  percent 
error  in  the  fit  to  the  strength  distribution  at  various  times. 

Figures  31  and  32  show  that  the  parameters  a  and  t  change  as 
the  power  n  of  the  decay  is  varied.  The  variation  in  t  is  always 
consistent,  increasing  with  increasing  n.  This  variation  is 
caused  by  the  requirement  for  the  distribution  to  reach  a  particu¬ 
lar  low  level  at  late  times.  If  the  assumed  decay  rate  is  faster 
(higher  n) ,  a  longer  delay  t  is  called  for  in  order  to  match 
the  experimental  data.  The  variation  in  a  with  n  shows  no  con¬ 
sistent  trends  and  seems  to  be  related  to  the  particular  set  of 
data;  similar  variations  in  a  are  seen  as  the  averaging  radius  is 
changed . 

Perhaps  the  most  interesting  parameter  to  examine  is  how  the 
rms  error  depends  upon  the  decay  power  n.  One  would  expect  the 
best  fit,  i.e.,  the  lowest  rms  error,  to  indicate  the  actual 
power  law  of  the  vortex  decay.  Figure  33  compares  the  power  n 
dependence  of  the  rms  error  for  all  ten  aircraft  types.  The  data 
show  a  consistent  trend  with  aircraft  size.  The  power  n  =  2 
gives  a  minimum  value  for  the  largest  aircraft  types  (B-747, 
L-1011,  and  DC-10).  The  B-707s  and  DC-8s  show  a  somewhat  slower 
decay  between  n  =  1.5  and  n  =  2.0.  The  B-727  data  indicate  n 
between  1.0  and  1.5,  while  the  B-737  and  DC- 9  data  show  no 
minimum.  The  values  of  n  for  the  largest  aircraft  are  likely  to 
be  the  most  accurate  representation  of  the  actual  vortex  decay; 
consequently,  the  value  n  =  2  will  be  used  in  subsequent  analysis. 
A  determination  of  the  value  of  n  depends  upon  making  measurements 
at  different  strength  levels  on  the  same  vortex.  The  most  accu¬ 
rate  decay  rate  measurements  occur  for  the  largest  aircraft  where 
the  initial  strength  may  be  four  times  the  MAVSS  detection  thres¬ 
hold  (see  Table  3).  The  smallest  aircraft  (B-737  and  DC-9)  may 
have  an  initial  strength  only  twice  the  detection  threshold.  The 
lack  of  a  well-defined  value  of  n  for  the  smallest  aircraft  is  a 
result  of  the  small  number  of  strength  bins  and  times  involved  in 
the  least-square  fit. 


FIGURE  33.  RMS  ERRORS  FOR  ALL  AIRCRAFT  TYPES 


The  data  on  least-square  fits  for  all  vortices  can  be  used 
to  assess  the  aircraft  dependence  of  vortex  decay.  In  order  to 
include  the  effects  of  both  t  and  a.  Equation  35  is  used  to 
calculate  the  times  tp  leading  to  the  same  hazard  probability 
F  =  0.2,  0.02,  and  0.002  used  in  Figure  30.  The  t  and  o  values 
for  n  =  2  and  10-meter  averaging  radius  were  selected.  Changing 
the  averaging  radius  has  little  effect  on  the  results.  The 
strength  threshold  is  set  equal  to  the  initial  strength  r^. 
Reducing  Fp  by  a  factor  of  two  would  increase  tp  by  /I  because  of 
the  n  =  2  decay.  The  dependence  of  tp  on  aircraft  wingspan  b  is 
shown  as  a  log-log  plot  in  Figure  34.  The  points  for  the  aircraft 
with  the  greatest  amount  of  data  lie  neatly  on  a  straight  line 
with  slope  .32  for  F  =  0.02.  The  points  lying  above  this  line 
are  for  the  three  aircraft  types  with  two  wing-mounted  engines 
and  the  B-707H  which  has  the  fewest  cases  of  all  the  types 
plotted.  The  presence  of  two  engines  on  the  wings  may  contribute 
to  abnormal  vortex  persistence. 

One  should  note  that  these  results  are  inconsistent  with  the 
conventional  way  of  normalizing  vortex  decay  to  the  time  bA/VC^ 
where  A  is  the  aspect  ratio,  V  is  the  airspeed,  and  CL  is  the  lift 
coefficient.  The  parameters  A  and  CL  vary  little  with  aircraft 
type.  The  landing  airspeed  V  may  be  at  most  20  percent  higher  for 
the  largest  jet  transports  compared  to  the  smallest.  Correcting 
for  V  would  increase  the  slope  in  Figure  34,  but  it  would  still  be 
less  than  half  the  unity  slope  expected  from  the  conventional 
normalization . 

In  principle,  the  parameters  of  the  least-square  fit  can  be 
used  to  examine  the  dependence  of  vortex  decay  upon  meteorological 
parameters.  In  practice,  the  use  of  two  free  parameters  t  and  a 
(n  =  2  is  fixed)  leads  to  a  poor  determination  of  each.  Figure 
35  shows  the  results  of  fitting  B-707  and  B-727  data  disaggregated 
by  vortex  number  and  wind  speed.  The  wind  speed  bins  were 
selected  to  have  approximately  equal  numbers  of  cases  (equal 
areas  in  the  x  and  y  wind  velocity  component  plane).  The  results 
show  that  the  parameter  t  generally  decreases  with  increasing 
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FIGURE  35.  LEAST-SQUARE  FITS  TO  DISAGGREGATED  B-707  AND  B-727  DATA 


total  wind  while  the  value  of  a  varies  haphazardly.  The  slower 

decay  of  vortex  2  is  sometimes  reflected  in  t  differences  and 

o 

other  times  reflected  in  a  differences.  The  way  in  which  the 

least-square  values  of  t  and  a  were  found  started  out  by  finding 

the  best  value  of  t  for  an  initial  value  a  =  10  seconds.  The 

o 

values  of  t  for  o  =  10  are  shown  in  Figure  36  for  two  wind 
speed  ranges  and  5  aircraft  types.  As  would  be  expected,  a 
single  parameter  fit  yields  values  that  are  more  consistent  in 
representing  the  differences  between  vortex  1  and  2  and  the 
effect  of  wind  on  vortex  decay. 

B . 3  SIMPLE  STOCHASTIC  MODEL 

The  least-square  fit  stochastic  model  of  Section  B.2  is 
inconvenient  to  use  for  predicting  hazard  decay  curves  since  it 
uses  a  measured  distribution  function  for  initial  strengths  and 
has  too  many  free  parameters  (t  and  a )  to  allow  an  unambiguous 
fit  to  the  data.  The  simple  stochastic  model  assumes  a  Gaussian 
distribution  of  initial  vortex  strengths  with  a  standard  devia¬ 
tion  equal  to  20  percent  of  the  mean  value  (see  Table  3) .  The 
two  parameters  tQ  and  a  are  reduced  to  one  by  assuming  tQ  =  3o 
which  is  in  the  middle  of  the  fitted  values.  The  value  of  n  is 
still  assumed  to  be  two.  Rather  than  conducting  a  complicated 
integral  over  the  initial  strength  distribution,  the  approxima¬ 
tion  is  made  that  the  hazard  probability  is  simply  reduced  by  the 
fraction  of  the  initial  strengths  lying  above  the  hazard  thres¬ 
hold.  Since  the  primary  purpose  of  the  model  is  to  extrapolate 
the  data  to  long  times,  this  approximation  should  introduce  no 
serious  errors.  The  resulting  hazard  probability  is  a  modifica¬ 
tion  of  Equation  33: 

F(I\j.,t)  =  \  [l  -  erf  ((Tj-r^)  / /I  ap)] 

x  £l  -  erf  ((t^T- 3a)  //7  a)j  (38) 
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Figure  37  shows  how  the  factor  describing  the  initial  strength 

distribution  affects  the  results  of  Figure  30.  The  straight  lines 

turn  downward  when  r4.  is  near  r*. 

l  o 

Figure  38  shows  one  hazard  probability  decay  curve  fitted  to 
the  form  of  Equation  38  by  a  steepest-descents  least-square  method 
which  minimized  the  deviation  in  F  (rather  than  InF  as  in  Section 
B.2).  An  optimum  value  for  a  was  found  for  each  of  the  six  r.j, 
curves;  the  a  values  are  listed  on  the  graph.  Only  measured  points 
for  F  <  0.5  were  used  in  the  fit  in  order  to  emphasize  the  decaying 
part  of  the  curves.  In  addition,  the  value  of  F  at  zero  time  was 
taken  from  the  data  rather  than  using  the  first  factor  in  Equa¬ 
tion  38,  since  the  resulting  fit  was  poor  for  the  cases  where  the 
calculated  initial  value  disagreed  significantly  with  the 
measured  value.  The  variation  in  a  shown  in  Figure  38  is  typical. 

The  values  are  roughly  constant  for  r|  between  the  detection 

2  1  2 
threshold  (about  50  m  /sec)  and  the  initial  value  (about  150  m  / 

sec).  Lower  values  are  obtained  for  r^.  outside  this  range.  One 

should  note  that  the  use  of  r^,  >  r  in  Equations  38  and  30  is  not 

strictly  consistent  with  the  original  model  since  the  vortex 

strength  was  supposed  to  be  constant  at  value  before  it  starts 

to  decay. 

The  results  of  fitting  the  decay  curves  in  Figure  38  were 
used  to  assign  a  single  value  of  o  to  describe  the  complete 
vortex  decay.  The  values  in  the  middle  of  the  constant  a  range 

7 

(e.g.,  r^.  =  75  and  100  m  /sec  in  Figure  38)  were  averaged. 

Since  the  purpose  of  the  selected  value  is  to  characterize  the 
low  probability  decay  region,  the  fitted  curves  were  adjusted 
slightly  to  give  a  better  fit  to  the  lower  part  of  the  decay 
curve.  This  change  generally  required  an  increase  in  a.  The 
final  value  selected  for  Figure  38  was  a  =  12.4.  Figure  39  shows 
how  the  calculated  curves  (using  both  factors  in  Equation  38) 
agree  with  the  measured  curves.  For  values  of  r,j,  at  or  below  the 
MAVSS  detection  threshold  the  measured  decay  is  more  rapid  than 
the  calculated  decay.  This  difference  could  be  the  result  of 
threshold  effects  on  the  measurements  or  it  could  mean  that 
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Note:  The  five  curves  represent  Tx  values  of  150,  100,  75, 
and  30  m^/sec  respectively  from  left  to  right.  The 
measured  data  are  plotted  with  thicker  lines  than  the 
fitted  curves,  the  numbers  below  the  curves  are  the 
values  of  a  for  each  curve. 


FIGURE  38.  FITTED  PROBABILITY  DECAY  CURVES  FOR  THE  B-707, 
10 -METER  AVERAGING  RADIUS,  AND  VORTEX  2 


FIGURE  39.  FITTED  PROBABILITY  DECAY  CURVES  USING 
TWO  PARAMETERS,  r'0  AND  a 


the  vortices  decay  more  rapidly  when  weakened.  In  any  case,  the 
calculated  values  are  a  conservative  estimate  of  the  duration  of 
low  strength  vortices. 

Figure  40  shows  the  values  of  a  selected  to  characterize 
the  vortex  hazard  decay  for  both  vortex  1  and  vortex  2  for  all 
ten  aircraft  types  and  four  averaging  radii.  The  value  of 
t  /a  =  3  used  in  the  model  gave  satisfactory  curve  shapes  for  all 
aircraft  except  the  DC- 9,  which  would  have  fit  better  with  a 
smaller  value,  and  the  B-747  and  L-1011  which  would  have  fit 
better  with  a  larger  value.  The  ratio  of  t  to  a  appears  to 
increase  with  aircraft  size.  In  some  ways  the  variations  in 
Figure  40  are  less  consistent  than  those  obtained  with  the  least- 
square  fits  of  Section  B.2.  Considerable  variation  with  averag¬ 
ing  radius  is  seen  here,  particularly  for  vortex  2,  while  little 
was  noted  before.  The  difference  may  be  due  to  the  fact  that, 
vortex  1  dominated  most  of  the  results  in  Section  B.2.  Figure  36 
also  indicates  more  drastic  variations  with  radius  for  vortex  2 
than  for  vortex  1.  The  general  trend  in  Figure  40  is  for  the 
smaller  radii  to  decay  more  slowly  (larger  a)  although  there  are 
obvious  exceptions  such  as  the  DC-8H.  The  differences  between 
the  DC-8  and  the  DC-8H  must  certainly  reflect  statistical  varia¬ 
tions  since  there  are  no  physical  differences  which  could  affect 
dependence  on  averaging  radius.  The  results  in  Figure  40  show 
too  much  variation  to  give  the  well  defined  size  dependence  of 
Figure  34.  However,  the  generally  slower  vortex  2  decay  for  the 
B-737,  B-707H,  DC-10,  and  L-1011  are  again  evident. 


B-26 


4 


0 


VORTEX  DECAY  PARAMETER  (Seconds) 


O  VORTEX  1 


AVERAGING  RADIUS  (m) 


FIGURE  40.  VALUES  OF  SIGMA  USED  TO  CHARACTERIZE 
THE  VORTEX  1  AND  2  HAZARD  DURATION 
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APPENDIX  C 
VORTEX  DECAY  PLOTS 


This  appendix  contains  the  vortex  decay  probability  plots 
for  the  B-747  (Figs.  41  to  48),  DC-10  (Tigs.  49  to  56),  L-1011 
(Figs.  57  to  64),  DC-8H  and  DC-8  (Figs.  65  to  72),  B-707H  and 
B- 707  (Figs.  73  to  80),  B-727  (Figs.  81  to  88),  B-737  (Figs.  89 
to  96),  and  DC-9  (Figs.  97  to  104).  Figures  65  through  80  are 
from  Volume  II  of  this  report  (Ref.  2) ,  but  are  repeated  here 
to  permit  comparisons  among  the  various  common  jet  transport 
aircraft.  Recall  that  vortex  1  is  the  first  vortex  to  be  detected 
by  the  MAVSS  sensors  and  is  the  downwind  vortex  (the  vortex  that 
might  translate  to  a  neighboring  parallel  runway,  for  instance). 
Vortex  2  is  the  second  vortex  to  be  detected  and  is  from  the  tip- 
wind  side  of  the  flight  path  (the  vortex  that  might  stall  near 
the  extended  runway  centerline).  The  plots  are  disaggregated 
into  vortices  1  and  2  and  winds  less  than  8  knots. 
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PROBABILITY  OF  DECAY,  B-747,  VORTEX  1  AND  2,  5-m  AVERAGING  RADIUS 
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FIGURE  43.  PROBABILITY  OF  DECAY,  B-747,  VORTEX  1  AND  2,  10-m  AVERAGING  RADIUS 


FIGURE  44.  PROBABILITY  OF  DECAY,  B-747,  WINDS  LESS  THAN  8  KNOTS,  10 -m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  B-747,  VORTEX  1  AND  2 ,  15-m  AVERAGING  RADIUS 


FIGURE  46.  PROBABILITY  OF  DECAY,  B-747,  WINDS  LESS  THAN  8  KNOTS,  15-m  AVERAGING  RADIUS 
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FIGURE  47.  PROBABILITY  OF  DECAY,  B-747,  VORTEX  1  AND  2,  20-m  AVERAGING  RADIUS 
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FIGURE  49.  PROBABILITY  OF  DECAY,  DC- 10,  VORTEX  '  AND  2,  5-m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  DC- 10,  WINDS  LESS  THAN  8  KNOTS,  5-m  AVERAGING  RADIUS 


FIGURE  51.  PROBABILITY  OF  DECAY,  DC- 10,  VORTEX  1  AND  2,  10-m  AVERAGING  RADIUS 


FIGURE  52.  PROBABILITY  OF  DECAY,  DC-10,  WINDS  LESS  THAN  10  KNOTS,  10-m  AVERAGING  RADIUS 


FIGURE  53.  PROBABILITY  OF  DECAY,  DC-10,  VORTEX  1  AND  2,  15-m  AVERAGING  RADIUS 


FIGURE  54.  PROBABILITY  OF  DECAY,  DC-10,  WINDS  LESS  THAN  8  KNOTS,  15-m  AVERAGING  RADIUS 


FIGURE  56.  PROBABILITY  OF  DECAY,  DC-10,  WINDS  LESS  THAN  8  KNOTS,  20-m  AVERAGING  RADIUS 


FIGURE  57.  PROBABILITY  OF  DECAY,  L-1011,  VORTEX  1  AND  2,  5 -m  AVERAGING  RADIUS 


FIGURE  58.  PROBABILITY  OF  DECAY,  L-1011,  WINDS  LESS  THAN  8  KNOTS,  5-m  AVERAGING  RADIUS 


FIGURE  59.  PROBABILITY  OF  DECAY,  L-1011,  VORTEX  1  AND  2,  10 -m  AVERAGING  RADIUS 


FIGURE  61.  PROBABILITY  OF  DECAY,  L-1011,  VORTEX  1  AND  2,  15 -m  AVERAGING  RADIUS 


FIGURE  62.  PROBABILITY  OF  DECAY,  L-1011,  WINDS  LESS  THAN  8  KNOTS,  15-m  AVERAGING  RADIUS 


FIGURE  63.  PROBABILITY  OF  DECAY,  L-I011,  VORTEX  1  AND  2,  20-m  AVERAGING  RADIUS 
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FIGURE  64.  PROBABILITY  OF  DECAY,  L-1011,  WINDS  LESS  THAN  8  KNOTS,  20 -m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  DC- 8  AND  DC-8H,  VORTEX  1  AND  2,  5-m  AVERAGING  RADIUS 


FIGURE  66.  PROBABILITY  OF  DECAY,  DC -8  AND  DC-8H,  WINDS  LESS  THAN  8  KNOTS,  5-m  AVERAGING 
RADIUS 
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FIGURE  68.  PROBABILITY  OF  DECAY,  DC- 8  AND  DC-8H,  WINDS  LESS  THAN  8  KNOTS,  10 -m  AVERAGING 

RADIUS 
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FIGURE  71.  PROBABILITY  OF  DECAY,  DC-8  AND  DC-8H,  VORTEX  1  AND  2,  20-m  AVERAGING  RADIUS 


FIGURE  72.  PROBABILITY  OF  DECAY,  DC-8  AND  DC-8H,  WINDS  LESS  THAN  8  KNOTS,  20- 
AVERAGING  RADIUS 


FIGURE  73.  PROBABILITY  OF  DECAY,  B-707  AND  B-707H,  VORTEX  1  AND  2 ,  5-m  AVERAGING  RADIUS 
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FIGURE  74.  PROBABILITY  OF  DECAY,  B-707  AND  B-707H,  WINDS  LESS  THAN  8  KNOTS 
AVERAGING  RADIUS 
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FIGURE  75.  PROBABILITY  OF  DECAY,  B-707  AND  B-707H,  VORTEX  1  AND  2,  10-m  AVERAGING  RADIUS 
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FIGURE  77.  PROBABILITY  OF  DECAY,  B-707  AND  B-707H,  VORTEX  1  AND  2,  15-m  AVERAGING  RADIUS 


FIGURE  80.  PROBABILITY  OF  DECAY,  B-707  AND  B-707H,  WINDS  LESS  THAN  8  KNOTS  20- 
AVERAGING  RADIUS 


FIGURE  81.  PROBABILITY  OF  DECAY,  B-727,  VORTEX  1  AND  2,  5-m  AVERAGING  RADIUS 


FIGURE  83.  PROBABILITY  OF  DECAY,  B-727,  VORTEX  1  AND  2,  10 -m  AVERAGING  RADIUS 


FIGURE  85.  PROBABILITY  OF  DECAY,  B-727,  VORTEX  1  AND  2,  15-m  AVERAGING  RADIUS 


FIGURE  87.  PROBABILITY  OF  DECAY,  B-727,  VORTEX  1  AND  2,  20-m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  B-727,  WINDS  LESS  THAN  8  KNOTS,  20 -m  AVERAGING  RADIUS 


FIGURE  89.  PROBABILITY  OF  DECAY,  B-737,  VORTEX  1  AND  2,  5 -m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  B-737,  VORTEX  1  AND  2,  10-m  AVERAGING  RADIUS 


FIGURE  92.  PROBABILITY  OF  DF.7AY ,  B-737  ,  WINDS  LESS  THAN  8  KNOTS,  10  -m  AVERAGING  RADIUS 


FIGURE  93.  PROBABILITY  OF  DECAY,  B-737,  VORTEX  1  AND  2,  I5-m  AVERAGING  RADIUS 
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FIGURE  94.  PROBABILITY  OF  DECAY,  B-737,  WINDS  LESS  THAN  8  KNOTS,  15-m  AVERAGING  RADIUS 
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FIGURE  96.  PROBABILITY  OF  DECAY,  B-737,  WINDS  LESS  THAN  8  KNOTS,  20-m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  DC' 9,  VORTEX  1  AND  2,  5-m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  DC-9,  VORTEX  1  AND  2 ,  10-m  AVERAGING  RADIUS 


PROBABILITY  OF  DECAY,  DC- 9,  VORTEX  1  AND  2,  15 -m  AVERAGING  RADIUS 
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GURE  104.  PROBABILITY  OF  DECAY,  DC- 9,  WINDS  LESS  THAN  8  KNOTS,  20 -m  AVERAGING  RADIUS 


APPENDIX  D 


DECAY  OF  AVERAGE  VORTEX  STRENGTH 


Figures  105  through  111  show  the  decay  of  average  vortex 
strength  for  the  DC-10,  L-1011,  DC-8H,  B-707H,  B-727,  B-737,  and 
DC-9,  respectively.  The  plots  for  the  B-747,  B-707,  and  DC-8  are 
in  the  main  text.  Figures  11  to  13,  respectively. 
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Averaging  Radius  =  5  m  (lower  line),  10  m  (middle  line),  and  20  m  (upper  line) 
FIGURE  105.  DECAY  OF  AVERAGE  VORTEX  STRENGTH  FOR  ALL  DC- 10  DETECTED  VORTICES 
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FIGURE  111.  DECAY  OF  AVERAGE  VORTEX  STRENGTH  FOR  ALL  DC- 9  DETECTED  VORTICES 
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